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PREFACE 

TniSBOOKis based upon a course of lectures given to engineers 
of the British Thomson-Houston Co. Ltd., and has been strongly 
influenced by their desires. Judging from them, it seems that 
most electrical engineers are aware that electrical transient 
problems can be solved by a mathematical method named after 
Oliver Heaviside; but the very fact that the method is mathe- 
matical puts it into a field in which many of them have practised 
little since their college days. The aim of this book is therefore to 
supply a link between the algebra and calculus of our youth and 
the circuit problems of our working hfe. 

The book has been written throughout with a bias towards the 
physical or engineering side; the question discussed is not ‘How 
do I solve this equation 1 ’ but ‘ How does this circuit work ? ’ 
There must indeed be a sound mathematical foundation, if the 
reader is to avoid the pitfalls which beset a too-practical ap- 
proach; but this necessary mathematics has been introduced 
gently, to reassure those readers whose practice in the subject is 
not of very recent date. As Mttle previous knowledge as possible 
has been assumed, but some knowledge which must be taken for 
granted has been summarized in the last four Appendices. 

The ground covered by the book is limited to the behaviour of 
linear circuits with lumped elements — a statement which is 
amplified and explained at the beginning of Chapter i. Within 
the field so limited the treatment is thorough, in the sense that 
I have not consciously omitted any important class of problem. 
Occasionally this policy has involved the discussion of somewhat 
abstruse ideas, such as the ‘stability’ considered in Chapter vi; 
especial care has been taken to present these parts in as simple 
a manner as possible. The Heaviside method has its own subtle 
difficulties, especially when it is applied to circuits which are not 
‘dead’ to start with. I have not always found these difficulties 
dealt with very clearly in the literature of the subject, so I have 
tried to ensure that the exposition of them is as simple and metho- 
dical as I could make it. 

No new step in the method has been introduced without being 
followed by one or more fully wq^ked-out examples of it; and all 
these examples have been dravm from real live engineering or 
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laboratory problems. The threads of the argument have been 
gathered together in a final chapter consisting entirely of worked- 
out examples, which between them illustrate most of the principal 
themes of the book. 

My tba.nks must first be accorded to my worldng colleagues, 
too numerous to name, who encouraged me to write the book and 
have helped me all through by their sustained interest and 
friendly criticism. The encouragement given by Mr E. H. Clough, 
C.B.E., Assistant Chief Engineer of the B.T.H. Co. Ltd., has been 
of especial value to me, and so has the help of Mr H. A. Price- 
Hughes, who has arranged for the proper drawing of the dia- 
grams. Among previous writers on this subject, my former 
teacher. Hr Harold Jeffreys, F.R.S., has helped me most, as any 
reader will discover if he reads Chapter i of his Operational 
Methods in Mathematical Physics alongside Chapters n and in 
of the present work. Finally, I am very grateful to my friend 
Mr J. W. Boag, who has checked all the algebra and arithmetic; 
to my wife, who has typed the whole manuscript; and to the 
staff of the University Press, for their care in the printing. 

G. W. C. 


January 1944 


CHAPTER I 


INTRODUCTORY 

1*1. The reader of this book is presumed to be already 
acquainted with electrical circuits; so he will realize that, in all 
their varied complexity, they are built up from a very small 
number of basic elements. Resistance, self-inductance, mutual 
inductance, and capacitance — these four exhaust the list. We 
shall be discussing circuits in which these four are variously com- 
bined; but the field which we shall cover will be narrowed in two 
ways. 

(1) The circuit elements will be linear. This means that the 
voltage across a resistance bears a constant ratio to the current; 
the voltage across an inductance bears a constant ratio to the rate 
of change of current; the voltage across a capacitance bears a 
constant ratio to the charge. In practice, this signifies the ex- 
clusion of elements like these: 

‘Non-linear^ resistance materials, such as are used in surge 
diverters; in these the current increases much faster than the 
voltage. 

Electric arcs. 

Inductances in which iron saturation is too great to be ignored. 

(2) The circuit elements will be lumped. Thus, we shall have to 
do with condensers, which may be treated as blocks of capaci- 
tance, or with coils, which may be treated as blocks of self- 
inductance; but not with such things as cables, every millimetre 
of which contains a tiny bit of self-inductance in the conductor, 
and also a tiny bit of capacitance between the conductor and the 
screen or sheath. A cable is a typical example of a circuit whose 
elements are distributed] the behaviour of such circuits constitutes 
a more advanced branch of our subject. 

It would be more exact to say ‘The circuit elements can be 
represented as being lumped’, for some circuits mth distributed 
elements are equivalent for most of our purposes to circuits with 
lumped elements. Thus, each turn of a coil has self-inductance, 
mutual inductance to the other turns, and resistance; yet the 
whole coil is equivalent to a block of setf-inductance in series with 
a block of resistance. It will be seen that the essential difference 
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2 INTRODIJCTOEY 

between sucli a coil and a cable is tbat, in tbe coil, tlie current 
is tbe same (at -a given moment) in every part; in a cable, it 
is not. 

When non-linear circuits and distributed circuits bave been 
excluded from our purview, the circuits which remain are neither 
uncommon nor contemptibly simple. In the field of power 
engineering alone, they comprise a great variety of power, testing 
and control circuits; exactly that class of circuit, in fact, whose 
behaviour under steady a.c. conditions is learned by every 
electrical engineering student. The reader will remember how, 
when each self-inductance i in a circuit has been treated as an 
impedance Ljo), and each capacitance C as an impedance l/OjV/;, 
the two components of both voltage and current in each part of 
the circuit may be calculated by a method just like the method 
used for n.c. circuits. We owe this method to C. P. Steinmetz 
(1865-1923), and by introducing it he enabled engineers to solve 
A.c. circuit problems without continual recourse to the drawing- 
board. 

This book treats of an analogous method for dealing with 
circuits under transient conditions, such as occur during the 
period of readjustment after a switch has been closed or opened. 
It is not much more difficult than the Steinmetz method, but it is 
by no means so widely known. It was devised by a man who used 
it on the problems of telegraphic transmission — Oliver Heaviside 
(1850-1925). He attained his results in a half-intuitive manner, 
which opened his methods to suspicion from the more orthodox 
mathematicians of his time. Heaviside’s obscure and dis- 
orderly writings are enlivened by sarcasm directed at these 
misguided rigorists, such as ' Whether good mathematicians, 
when they die, go to Cambridge, I do not know’. How^'ever, 
Cambridge eventually returned good for evil by the hand of 
T. J. I’A. Bromwich (1875-1931), who showed that Heaviside’s 
results reaUy could be soundly based on mathematical reasoning. 
The same thing was done simultaneously by K. W. Wagner of 
Berlin; independently, for the date was 1916, and the native 
countries of the two workers were at war. 

1*2. A pair of simple examples will show the parallelism be- 
tween the Steinmetz and the Heaviside method. 

(1) An alternating voltage F is applied to the circuit shown in 
Pig. 1. What wdU be the voltage v across the condenser? 
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Representing the condenser as an impedance 1 ICj(t>j we see that 

V ___ IjCjo) 

V E+l/Cjo) 

- 1 

1 -{-jMCo) 

l-jMCo) 

“ l + (i2Ca>)2’ 

so V has the component F/1 + [RCiof in phase with F, and the 
component ~ RO(i)Vjl + {RC(j)Y in quadrature with it. 



Fig. 1. Fig. 2. 


(2) In the circuit of Fig. 2, the switch S is closed, the condenser 
having initially no charge. A d.c. voltage F is thereby applied to 
the circuit. Find how the voltage v on the condenser increases. 

In the Heaviside method of solving this problem, the con- 
denser is regarded as an impedance l/(7p. Then 

V _ IjCp 
f^R+ljOp 

1 

RCp + 1 

= 1 L. 

f + llRC 

'P 

In other words, v =V 

p+llBC 

Reference to a table of standard forms, such as is given in 
Appendix I, shows that 

V 

p + ljEC 

is to be interpreted as meaning 

F X 
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where t denotes time from the instant of closing the switch. 
Therefore 

V = F{1 — 

which gives the value of v at every instant after the time t = 0 
when the switch was closed. 

Simple though this process is, it is unwise to learn it by rote, 
as a child learns long division. The foundations of it need to be 
understood, if it is to he applied without fear of pitfalls. To ex- 
plain these foundations as simply as possible is the task of the 
next few chapters; then we shall go on to some practical applica- 
tions. 



CHAPTER II 


TRANSIENT CONDITIONS IN THE 
SIMPLEST CIRCUITS 

2T. The Differential Equation of a Circuit. 

Consider further the circuit of Fig. 2. At every instant of time 
(denoted by t), after the switch has been closed, the charge q on 
the condenser is related to its voltage v by the equation 

q = Gv. 

Therefore the rate of increase of o' = C x the rate of increase of v, 
since (as has been stated once for all in § M) G does not vary: 
in other words 

dq _ dv 
dt dt ' 

But the rate of growth of the charge is equal to the current which 
is flowmg into one plate of the condenser and out of the other. 
Passing through the resistance i?, this current gives rise to a 
voltage Bdqjdt, or RGdvjdt. Now, at every instant. 

Resistance voltage + Condenser voltage = Apphed voltage F. 

Therefore BG ^ + v = V. 

dt 

This equation, which expresses the behaviour of the circuit, is 
called the differential equation of the circuit; the solution oi this 
equation is a relation between v and t, which enables the value 
of V at each instant t to be calculated. If the condenser started off 
already charged to a voltage v^, we should expect to find a dif- 
ferent solution from that corresponding to an initial state of 
complete discharge; and in fact no problem whose differential 
equation involves derivatives of the first order (dvjdt, but not 
d^v/dt^, dHjdt^, ...) can be solved without a knowledge of one such 
starting-point (or its equivalent). 

2-2. Solution of the Differential Equation. 

Before proceeding to solve the equation 

■n„dv 

BG-^ + v = V, 
dt 
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let ns rewrite it in a form which can embrace many other problems 
besides this one. To begin with, the applied voltage V was taken 
as a D.c. voltage for simplicity; but it might have been an alter- 
nating voltage, or even some more complicated shape. All these 
are covered if, instead of the constant V, we write a function of t, 
say f(t). Again, we might have wanted to calculate charge, or 
current, instead of voltage; we cover all these possibilities if we 
write a non-committal u to stand for any of them. So the final 
form of the equation will be 


a~ + bn=m. 


where a and b are constants. We will seek that solution which 
makes u = when t~0. 

We shall use the letter Q to denote the operation of integrating 
with respect to t, between the limits 0 and that is, Q<p(t) is a 

shorthand way of writing f (p{t)dt. Therefore, for instance, 

Jo 

Qt = t^l2, 

Q^i = 

and so forth. We perform the operation Q on every term of our 
equation, so as to get 

+ = C/(<), 

and proceed to consider the meaning of each of these terms. 



= a{u — UQ). 


(2) Q{f>u) = f budt 

Jo 

= bQu. 
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We do not know as a function of so we leave this term 
in. the last-named form. Note that it is permissible to change the 
relative positions of Q and the constant b. 


Qm = f 

Jo 


which can be evaluated if necessary, but will be left for the 
moment in its original form. 

Rewriting the equation in the light of these statements, we 


obtain 


a{u--iiQ)-TbQu = Qf{t), 


U = na+-Qfit)~^Qu. 

Ob 


We shall now start to VTite a instead of bja. 

The first two terms of the right-hand side are known, but not 
the last term. However, for u in the last term we can substitute 
the expression 

'i^Q + ~Qf{t)-ccQu, 

Of 

so as to get 

M = ^ Qf{t) - aQ {mo + ^ Qf{t) - ccQu]. 

d d 


If the same substitution is performed repeatedly, we get 

u = Uq + ^ Qfit) - ccQ + - Qf{t) -aQ[uo+^Qf(t)- ... I 

so that the unknown term is, as it were, pushed farther and farther 
away. 

Now aQaQu^ is the same thing as where means that 

the operation Q must be performed twice. Therefore, collecting 
together all the Uq terms and all the/(^) terms, we get 

% = {l-aQ-i-a^Q^-a^Q^+ ...)uQ+{Q-ccQ^ + a^Q^- 

Every term of these expressions can be evaluated by integrating 
once, twice, thrice, four times, so here is an equation which 

makes it possible to work out in terms of 
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The bracket leads readily to a simple result. For 

qcQuq = a f u^dt = UQ.ai, 

Jo 

a^Q^UQ = u^atit ~ 
a^Q^UQ = Uq.ocH^IZI, 

and so on. Therefore 

(Qr2f2 ^ 3#3 \ 

+ — 31’^*“)* 

We recognize this series as the exponential one, whose sum is 
Therefore 

u = UQe~^ + {Q-aQ'^ + a^Q^- ...) j^“/(0 • 

There is no obvious way of finding the sum of the second series 
without knowing the exact form of f{t)] so at the present stage we 
can only hope that it will usually prove as recognizable as the 
first one. Mathematicians might be seized at this point with a 
horrible misgiving that the series may be ' divergent ’ (or meaning- 
less); but we can be fortified by the reflection that our problem, 
being a physical problem, must have a solution which means 
something. 


2-3. Example. Solution of a Circuit by the Q method. 

Let us apply this method to the circuit of Fig. 2, whose differential 
equation was shown to be 

but let us suppose that, prior to the closing of the switch, the 
condenser was charged to a voltage Vq . 

The equation may be written 


dv 

dt 


+ av — aV, 


where a = IjEC. The operation Q, performed on every term, 
gives 

v — VQ + aQv = (xQVy 


V = VQ-\-aQV — aQv; 


or 
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and the method, already explained leads on to 

V = Vo + aQV-ocQ{v^ + aQV-ciQ[Vo + ccQV -...]} 

= il-ccQ + oc^Q^-...)vo + {o!,Q-a^Q^ + a^Q^-...)V 
= Voe~=‘i+(ciQ-a^Q^ + a.^Q^-...)V. 

We observe with, joy that the second term also leads to a recog- 
nizable series, since F is constant ; namely 


Vioct- 


ocH^ 
‘ 2 ! 


aH^ 


-t- 


s.'-). 

which is an exponential series shorn of its first term, and equals 
So the solution comes to be 

where ci — IjRG, This shows that the 
relation between v and t is of the form 
shown in Fig. 3; starting at the initial 
value ^’Q, the voltage climbs up to the 
full value F. 

I hasten to assure the reader that the Heaviside method does 
not involve this ponderous process, except in the earliest stage of 
learning it. 



2*4. Behaviour of the Operator Q. 

The method just described is called an operational method, 
because it is carried out by means of the ‘operator’ Q. (There is 
nothing more abstruse about the use of the symbol Q to denote 
the operation of integrating, than about the use of the symbol x 
to denote the operation of multiplying.) Hitherto we have treated 
Q with respect, remembering always that it is not a number but 
an operator; but now we shall discover that great abbreviation is 
possible if we appreciate that it can in fact be treated very nearly 
as a number. 

In solving the equation 

we reached almost immediately the stage 

u = Uo-\--Qf{t)-aQu. 

a 
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Suppose we had solved this equation for m, treating (3 as a number. 
We should have obtained 

Now if Q were a number, 1/1 + aQ could be expanded by the 
Binomial Theorem (see Appendix II), with the result 

= l-ocQ + a^Q^-oc^Q^+ .... 

l+aQ 

This would have led us immediately to the equation 
u = {l -aQ + a^Q^- ...)uQ + (Q-aQ^ + oc^Q^- ...) 


which we got before, with much more labour. The method leads 
to the right answer, so there must be something in it. Why are we 
justified in treating Q in this way? and do any pitfalls lurk to 
entrap us when we do it? 

If u and V are two functions of t, 

Q{u + v) — J {u + v)dt 

== f udt+ f vdt 
Jo Jo 

= Qu^Qv, 

just as if Q were a number. Similarly 

QmQn^ = Q^+nu^ 

and, if a is a constant, 

Q{au) = aQu. 

But Q{uv) is not equal to uQv: for instance, Q{t.t) = Qt^ = f jZ, 
whevesbSt.Qt = t^/2. 

Thus we may say that Q obeys the ordinary laws of algebra, 
except that it must not change places with a variable. We may treat 
Q as a number, provided that we take care that its position 
relatively to variable terms is not altered. For example, w^e must 
not change Q{t sin t) into t . Q sin t, for they do not mean the same 
thing. 
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There is therefore a Justification for the rapid transition from 

n = + ~ Qf{t) - aQu, 

(h 


or 


(l + aQ)u — Uq + — Qf{t), 




Neither in the division of both sides by (1 + a$), nor in the re- 
placement of IJl+aQ by its binomial expansion, is there any 
violation of the rule that Q must not change places with a variable. 

2*5. Heaviside’s Operator p. 

We are on the verge of being able to write down a simple rule for 
solving any circuit of this type; but before doing so we have to bid 
farewell to our operator Q, since it is not in general use. The 
operator used by Heaviside and his successors is denoted by the 
letter p, and we shall define it by the equation 


P = 1/Q. 


For the moment, the reader may think of this change as being 
simply a change of notation; for reasons not yet revealed, we 

prefer 1 jp to Q, as a notation for the operation So far, we 

Jo 

have attached a meaning to Q, <3^ ..., but not to IjQ, 
therefore we can attach a meaning to 1/^, 1 jp^, . . . , but not (as yet) 
to p, p ‘^, .... Q is more convenient in proving the method, but p 
is more convenient in applying it. 

All that has been said about the legitimacy of treating Q as a 
number applies equally to p ; it obeys the ordinary laws of algebra, 
except that it must not change places with a variable. 

2*6. Rules for Solution of a Simple Electrical Circuit. 

The formal solution 

now takes the form 

u = Uq + — h fl /(*)"] 

p + a, “ P + J 


ap 

ap + b 


ap+b 


m- 
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We may therefore write down the following rules for the solution 
of any circuit which is simple enough to lead to a differential 
equation of the form 

(It is not necessary to define what circuits lead to equations of 
this type, since we shall almost immediately be proceeding to 
more complicated circuits.) 

(1) Write down the differential equation of the circuit for the 
quantity desired (e.g. voltage), in the form 

6m =/(«)• 

(2) On the left-hand side, replace djdt by p, 

(3) To the right-hand side, add the result of dropping the 
-z^-term from the left-hand side and changing {pu) into {pu^). 

(4) Solve the resulting equation algebraically for u (treating p 
as a number), and evaluate the result by expanding in powers of 

Ijp and interpreting 1/p as f 

Jo 

The first stage may afford difficulty to some. Later on we shall 
show that stages (1) and (2) can both be replaced by a process 
closely resembling the Steinmetz method for solving a.c. circuits 
under steady conditions. For the present, we merely recall that 
Voltage and current in a resistance R are related by 

V — Ri\ 

Voltage and current in an inductance L are related by 
V — Ldijdt] 

Voltage and current in a capacitance G are related by 
i = Odvjdt; 

while in a mutual inductance M, whose primary and secondary 
coils have self-inductances ii, 1 / 2 , 

Primary voltage 

dt dt 

and Secondary voltage — + — . 

dt “at ' 
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These equations of course involve certain assumptions about the 
relative directions of current and voltage — ^if one is reversed, a 
negative sign will creep in. 

2-7. Standard Operational Forms. 

We do not in practice carry out every time the expansion in 
powers of 1/j). We do it once for all, and record the result for 
future reference as a ‘Standard Form’. (See Appendix I.) It will 
be seen that the solution of this type of circuit involves two such 
forms; namely 

P 

— — operating on a constant, Un, 

jJ ~T (Z ^ 

, 1 

and — ; — operating on a function of t. 

p i" a 


We proceed to find out the general meaning of these. Now 


P 


.1 = 


..1 


+ l + ciQ' 


-c&i 


as has been seen already. Also 

equation a 

du 


f(t) is that solution of the 


dt 


-\-au = f{t) 


which makes = 0 when ^ = 0. To find this, multiply both sides 
of the equation by we get 


or 





d 

dt 




or 

since % = 0 when ^ — 0; 



or 
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Thus -^/(O = 

p + oc Jo 

or, in words, 

'The operation l/_2> + a means: Multiply by integrate the 
product from 0 to t, and multiply the result by 

This last result gives a simple expression for something which 
hitherto we have known only as an infinite series, namely 

{Q - olQ^ + _ . . .)/(^). (See § 2-2.) 

Instead of having to integrate /(i) once, twice, thrice, ..., hoping 
to discover some rhyme or reason about the process, we now have 
to perform a single integration only. In Appendix I will be found 
the result of performing the operation 1/^ + a on certain com- 
monly encountered functions of t. 

This chapter will be concluded with some practical examples. 

2*8. Examples. Solution of the Simplest Type of Circuit. 

Example 1. The field winding of an alternator has a self- 
inductance L henries, and a resistance R ohms. It is passing a 
current of amps, from a n.c. supply (Fig. 4). By the operation 
of a field suppression circuit, the winding is suddenly discon- 
nected from the supply and short-circuited on itself. Find how 
the current now varies. 

The operation is carried out by moving the switch shown in 
Fig. 4 from position A to position B; the circuit must be closed at 
B before it is opened at A. There will also be a parallel circuit 
(not shown) for receiving the current which is still wanting to 
flow from the d.c. supply. 

At any instant t seconds after the switching, suppose the 
current to be i amps. The inductive voltage in the field winding is 
Ldijdt; the resistance voltage, Ri, Since the winding is short- 
circuited, 

Voltage on i-f voltage on i? = 0, 

that is, L^ + Bi = 0; 

at 

and, when ^ = 0, i = i^. 

This is the differential equation of the circuit. The p-equation, 
by the rules given in § 2*6, is 


{Lp + R)i — Lpi^. 
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_ Lp 

Zrp + i?*' 


= . Li 

Thus the current collapses according to an exponential law. It has 
fallen to half its original value when t has the value given by the 

equation 

p~miL _ 1 

— 2J 

or = 2, 

or RtjL = loge2 = 0-693, 

or t = 0-69ZLIR. 


SWITCHING RESISTANCE 

(opened when contact B 

13 CLOSED) 





B o 

o 

D.C 

Supply 



Fig. i. 



Pig. 5. 


Example 2. A transformer, whose windings have self-in- 
ductances henries and a mutual inductance M henries, is 

connected on its secondary side to a resistance load of R ohms; 
the resistances of the windings are neglected (Fig. 5) . A switch is 
closed on the primary side, connecting the transformer to an a.c. 
voltage supply of peak value V and frequency i3/27r, at the 
instant of voltage maximum. Before the switch was closed, 
the transformer was dead. Find the voltage across the load 
resistance i2. 

If the primary and secondary currents are % and amps., we 
have 


y. dt-t 


: Voltage applied to primary terminals 


= V cos Qt, 
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since tlie voltage has to be a maximum when t 
sidering the secondary circuit, 


since there is no source of voltage in this circuit (apart from the 
transformer voltage, which is represented by the Mdii^jdt term). 
We are only interested in the secondary circuit; so we eliminate 
by multiplying the former equation by Jf/Li and subtracting 
it from the latter. This gives 


Replacing by vjR, where v is the desired secondary voltage, we 
obtain 

/L^L^~M^\dv 


This is the differential equation. 
The ^-equation is therefore 


for, with ^ = 0 initially, there is no Vq term to add. Therefore 


aMV' 


where 


Reference to Appendix I (form 23) shows that this is to be 
interpreted as 

MV a , 

The behaviour of in a typical case is shown in Mg. 6. The last 
term in the expression is a transient one, and when it has died out 
we are left with an alternating voltage only. 
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Tlie constant a may be written RjL^ where 




w 


L is the short-circuit reactance of the transformer, viewed from, 
the secondary side; and ( — MVjLi) peak value of the e.m.f. 
referred to the secondary side. 



If the steady a.c. voltage across R is worked out by the 
Steinmetz method, assuming series inductance L and e.m.f. 
( — MV jLj)y the result will be found to be 


MV a 


(a-jQ), 


which is of course equivalent to the first two terms of the expres- 
sion which we have worked out for v. 


CET 
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CHAPTER III 


TRANSIENT CONDITIONS IN MORE 
COMPLICATED CIRCUITS 

31. Significance of the Operator p. 

Before going on to more complicated ckcnits, we pause to ask a 

question: ‘If Q = J dt, and p == I jQ, what is p? What physical 

meaning, if any, does it possess? ’ 

Since p replaces djdt in the equations, it is tempting to suppose 
that p means djdt; and indeed that meaning is often attributed 
to it. But this interpretation may easily lead us into trouble. For 
instance, in Example 1 at the end of the last chapter (§2*8), the 
auxiliary p-equation was 

(ip + i?)i = 

and if p means didt, the right-hand side is zero, so that we are led 
to the odd conclusion that the decay of current in a field sup- 
pression circuit is not affected by the current you start with ! 

A truer insight may be obtained, if we observe that the equation 

P = 1IQ 

means that p is the operator which undoes the operation Q. Does 
djdt do thm'i Assuredly 

-/«); 

but = jf'{t)dt 

=/(«)-/(0). 

Thus djdt undoes Q, if, and only if, it acts after Q (apart from the 
special case/(0) = 0). In other words, p can only be identified 
with djdt if the differentiations are performed after the integra- 
tions. The wisest course, however, is to remember that p means 
l/Q, and forget its identification with djdt. Only seldom is the 
latter necessary. - 
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3-2. Circuits of the Next Stage of Complexity. 

In Chapter n we considered exclusively such circuits as would 
lead to differential equations' of form 

du , 

111 practice, this meant simple combinations of resistance with 
either inductance or capacitance. If we consider two simple 
resistance-inductance circuits with mutual inductance between 
them, we have a form which makes a convenient stepping-stone 
to a general theory. 

In the circuit of Fig. 7 the switch is closed at time ^ = 0 so that 
the circuit (hitherto dead) is connected to a I -volt supply. How 
do the primary and secondary currents behave? 



Fig. 7. 


Let the currents be and then, just as in Example 2 at the 
end of the last chapter (§2*8), the differential equations of the 
circuit are 


_ d'l-t d'ba , 

du 

-^ + 3 ^ + 2 *, = 0 . 

Now if the initial values of and were and the auxiliary 
:p-equations (written down by rules exactly analogous to those 
laid down for a single variable) would be 

and -f ( 3 ^ -}- 2) -{- Spigo- 

In this particular case, ^20 zero 5 so the p-equations are 

simply'' 

(2p-f- i)ii+pi2 = 1, 

pii^ -f (3p -}- 2) ^2 = 0, " 



2-2 
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Solving for treating p as a number, we get 


-* 3j9-f 2 

(2j) + 1) (3j? + 2)—p^^ 

^ Zp + 2 ^ 

5p^ + 7^ + 2 * 

As it stands, this expression is not a standard form. The way to 
interpret it is to split it into ' partial fractions ’ which are standard 



forms. Rules for this will be given presently, so we will simply 
write down the result here. The above expression for is the 
same as 

. _ 3p + 2 

(5p + 2){p+l) ’ 

and the reader may verify that this is algebraically equivalent to 

i = ^ V \ 1 

^ I 3p + (f) 3p + lj-^ 

= 1 — I e-2^/5 1 

Similarly ^2 = 

The currents therefore vary in the manner shown in Fig. 8. As 
is to be expected, tends towards the d.c. value of 1 amp., 
whereas settles down to zero when the effects of the switching 
operation have died away. 

3*3. The Differential Equation of Any Linear Circuit. 

In solving this last example, we wrote down a pair of differential 
equations which contained both and ^ 2 ? and proceeded immedi- 
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ately to solve them by the p method. But we might instead have 
got rid (for instance) of ^ 2 , and so obtained a differential equation 
for % only; this would, in fact, have been 




We observe that this is an equation of the general form 


a 


d^u du 


+CU =f{ty. 


and we recall that the circuits of the last chapter led to the form 

du , , 

=/(<)• 


It begins to dawn upon us that the differential equations associ- 
ated with any circuit with linear lumped elements are likely 
to be of the form 


d^u ,d'^~'^u ^du . 


and indeed this does prove to be the case. 

These equations are all of the form known as linear differential 
equations with constant coefficients; and if the highest-order 
derivative occurring in the equation is d^ujdt^, the equation is 
said to be of the ?^th order. Our chief interest is in the circuit, not 
in the differential equation; but it saves a lot of circumlocution 
if we write ^ a circuit of the wth order instead of ' a circuit whose 
associated differential equations are of the form 


d^'U ^ du „ . , , 

u-T7r + ...+Z — =f{t)\ 


dP^ 


dt 


It now becomes plain that the circuits solved in Chapter ii were 
circuits of the first order; and that the circuit of Fig. 7 is a second- 
order circuit. 

3-4. Further Standard Forms. 

These higher-order circuits are capable of leading top-expres- 
sions of other forms than the two already discussed: as, for 
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instance. 


V 




. 1. These can all be interpreted by expansion in 


powers of ljp\ for instance, 


P 


.1 = 




.1 


/_ a? 

1 1 — 5 + — i — 
\ p^ p^ 

aH^ 



= cosat. 


A list of such standard forms will be found in Appendix I. Some 
of the more complex ones may be neatly derived from simpler 
ones like this: 

We know, for instance, that . 1 = 

p-i-a 

Differentiate both sides of this equation with respect to a, 
treating p and t as constant. We get 


1 = 

{p + a,f 




or 


P 

{p + a)^ 


. 1 = te’^, a new standard form. 


It may seem strange to treat the one constant in the equation as 
a variable; however, the method leads to the right answer, so it is 
not so unsound as it looks. 


3-5. Method of dealing with Operations upon_/’(t). 

It will be seen that most of the standard forms in the list are con- 
cerned with interpreting an expression of the form F(p) . 1 ; but 
we know from experience with first-order circuits that we may 
also have to interpret expressions like F{p).f{t). In very many 
circuit problems, however, /(^) is an elementary function such as 
a sine wave. For such functions, an operational equivalent may be 
immediately picked out from Appendix I; so that, if the equi- 
valent of 

fit) is f>{p).l, 
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we may immediately write 

and the combined function F{p)(j){f) can be split up into a number 
of standard forms of the kind in which a function of p operates 

upon 1. 

For instance, if it were necessary to interpret 


3 1 


• sin 


we should use Appendix I to write sin 2t in the form 

2p 


P^+4: 


. 1 . 


Thus 


1 


2p^+l 


sin 2t = 


Zp 


{p^+ l)(p2 + 4) 


1 


= / ^ P \ 
= sin ^ — I sin 2t. 


3-6. Heaviside’s Shifting Theorem. 

A theorem which is sometimes used is here mentioned for 
reference only, in case the reader should find it named in the 
literature. If a function of ^ operates upon a function of t having 
an exponential factor, it is possible to remove the exponential by 
the use of the identity 

F{p).{e‘^^f(t)} = e‘‘*F{p + a).{fit}}. 

This can sometimes lead to a short cut, or to a method of deriving 
a fresh ^standard form’. But if the reader is prepared to take 
'standard forms’ for granted, he can safely leave it unlearnt. 


3-7. Formation of Partial Fractions. 

The reader will have gathered that it is necessary to split up 
complex algebraic functions of p into a number of expressions 
each simple enough to be a standard form. The following rules 
will enable this to be done. 
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We will suppose that we have to interpret such an expression as 

+...+ 61 ^ 9+60 ^ 

4- . . . + a# + Uo ‘ 

We always find in circuit- problems that m is not greater than n. 

It is shown in Appendix V that algebraic expressions, such as 
the top or bottom lines of this fraction, can be split into real 
factors: 

either single linear factors {p + a), 

or repeated linear factors {p + , 

or single quadratic factors 4 ap + 6 ), 

or repeated quadratic factors (p^ 4 cp 4 df. 

In the last two cases < 46, < 4d: for otherwise the expressions 

could be further split up into real factors (p 4 7 ) (p 4 . Quadratic 
factors are only introduced when the linear ones are not real. 
The partial fractions are found as follows. 

(1) Rewrite the bottom line of the operational fraction in its 
real factors, thus: 

4 4 . . . 4 61 P 4 60 

ajp4a) ... (p4y5)’*... (p^4ap46) ... (p^4cp4d)«... ’ 

(2) Equate this expression to the sum of a string of fractions: 

for each single linear factor - assume a fraction of 

Ajy {p + a)... 


form 


p4a 

for each repeated linear factor 


of form 


Bp Cp 


4 


;4 ... 4 ; 


Op 


ip+j^r 


assume fractions 


p+fi {p+pf (p+^y 

for each single quadratic factor 

fraction of form 

p^+ap+b 

for each repeated quadratic factor 
fractions of form 

Lp^ + Mp ^ Np^ + Pp 


ip^ + ap + b) 


assume a 


i + 


p^ + cp + d (p^+cp + df 


4 ... 4 


(p^4cp4c^)®.. 

/Sp^4jPp 


assume 


(p^4cp4e^)’^ 
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and in every case assume also a constant term JT. Then it will be 
found that it is possible to choose numbers for A,B,G, 
so that the sum of the string of fractions is equal to the original 
expression 

+ ■ ■ ■ + 6ij) + fep ^ 

®«,(i5 + a) ... (25 + /?)’' ... + + ... {’p^ + cp + dy...' 

There will be just enough disposable numbers -4, S, ... for this; 
neither too many nor too few. 

There exists a general formula, known as Heaviside’s Ex- 
pansion Theorem, for writing down the function of t which 
emerges when any algebraic function of p operates upon I. In 
practice it is less convenient than the method just given, and is 
only mentioned here in case the reader should find it named in the 
literature and desire to know what it is. 

3-8. Examples of the Formation of Partial Fractions. 
Example 1. The circuit of Fig. 7 led us to the expression 

i = 1 

^ 5p^ + lp + 2' 

Hp+i)(p+i): ■ 

Let - - s A + + 0. 

5(f> + f)(2’+l) P + ^s 2> + l 

Then (multiplying both sides by 5{p + l){p + 1)), 

3p + 2 = QAp{p + 1 ) + 5Bp{p + 1) + 5G{p + 1) (p + 1 ). 

The simplest way to get A and B, in such a case as this where the 
factors are linear, is to give p particular values, which cause the 
disappearance of every group except one. Thus 

p = - 1 gives -f + 2 = 5J:( - 1) (I), whence 4 = - 1; 
p = — I gives - 3 + 2 = 5JB( - 1) ( — 1), whence B = — 

The constant term C can always be obtained by writing p = 0; 
here, this gives (7=1. Therefore 

^~[~3p + % 3p + l^ /■ ’ 

which is interpreted in the way already indicated. 
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Example 2. In § 3-5, we had to interpret the expression 


Following the rules, we assume 

(i52+l)(35^4-T)^ x>^+l ^ i>2 + 4 

Writing = 0 shows that jB = 0. For the rest, multiply up, so as 
to obtain 

3jp = 4- Bp) (p2 + 4) -f ((7^2 ^ ])p'^ (1^2 

The right-hand side contains terms in and p. Equating 

the coejfjficients of these with the corresponding coefficients on the 
left, we obtain 

A + G = (}, 

B + D == 0, 

4cA -j- (7 = 0, 

45-fD = 3. 


The first and third give ^4 = (7 = 0; the second and fourth give 
£ = 1, D = - 1, Therefore 

(p’ + W + i)-' ■ '■ “ ’''‘■O"- 

3-9. Circuits of the Second Order. 

We have seen that there is a class of circuit whose differential 
equations have the form 


W M/ , U/M/ ^ 

Just as, in the first-order circuit, we needed to know the starting 
value of one quantity, so in the second-order equation we must 
know two starting values (or their equivalent). For instance, it is 
sufficient if we know that u = Uq, dnjdt = u-^, when t = 0. 

The equation is easily split into two of the first order, by taking 
a new variable v = dujdt. Then 


dv , 

a-j^+bv + cu =f(t), 


du 

dt 


'V = 0, 
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The auxiliary j^-equations, written down in the usual way, are 

{ap + 5) V + = f{t) + 

pu — v= pUq. 

Eliniinatiiig v algebraically, we find 


(ap^ + bp-{-c)u = f{t) + {ap^ + bp) Uq + apuj^. 


This can be interpreted by splitting each term into suitable 
partial fractions. We shall not go through the process here, for it 
will be performed in several examples at the end of this chapter. 

3-10. Circuits of Any Order. 


In exactly the same way we can solve a circuit of order n, given 
that, when t = 0, u = Uq, dujdt ~ d^ujdt^ = etc. It soon 
appears that the p -equation may be written down from the dif- 
ferential equation by the following rules. 

( 1 ) On the left-hand side, replace djdt by p. 

(2) On the right-hand side, leave/(^) unchanged, but add terms 
in W0, u-i, 'leg, . . . obtained in the following manner: 

Uq term. Drop the term which does not contain p from the 
left-hand side, and write Uq for u\ 

term. Divide the last named by p, drop the term which 
now does not contain p, and write for 

u^ term. Divide the last named by p, drop the term which 
now does not contain p, and write for up, 

and so on. 

Thus, the third-order differential equation 


d^u 


+ b 


d^u 


du , 
+ c^ + du 
at 


m 


leads to the p -equation 
{ap^’Php^‘\-cp-{-d)u 

= f{t)-^{ap^pbp^-pcp)uQP{ap^-^bp)u^^apu^, 

The extension to simultaneous equations (involving two or more 
dependent variables u, v, w , involves no new principle but 
is : rarely required.' 
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3-11, Examples. Circuits of Order Higher than the First. 
Example 1. This circuit resembles that ^ 
shown in Fig. 2, except that inductance L 
has been added (Fig. 9). The closing of the 
switch 8 applies to the circuit aD.c. voltage F ; 
how does the condenser voltage v behave, as- 
suming that it was zero to start with? 

We have seen (§ 2*6) that the current and 
voltage in a condenser are related by the 
equation 


Fig. 9. 


. ^dv 
% — G — . 
dt 


di 


Therefore the inductance voltage = = LG 

ttb 


the resistance voltage 
the condenser voltage 
and the diiBferential equation is 


Ri = EC 


= v: 


dv ^ 

Tf 


rr^d^V jy^dv -rr 


We must now consider what are the initial values of v and dv/dt. 
The condenser is initially discharged; its voltage cannot give a 
leap when the switch is closed, for that would involve passing an 
infinite current through L and B; so Vq = 0. Again, dvjdt is pro- 
portional to the current, and this cannot change suddenly because 
the voltage across L would then be infinite; so Wj, the initial value 
of dvjdt, is zero also. Hence the ^-equation is 


and 


{LCp^^RCp+l)v= F, 


V — 


2 , -S 

'P^ + -yP + 


LG 


V_ 

LG' 


The solution of this may take one of two forms, according as 
B^ is greater or less than 4^/0 (see Appendix V). If B^ > ^LjG, 
the denominator has real factors (p -i- a) (p •+■ §). If B? < iLjG, the 
factors are complex, and the denominator may be written in the 

form (pH- 7)24- 
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= ^ V 

(j5 + a)(p+yff) ’ 
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for — IjLO. Assuming partial fractions 

find 


Ap 

p + a 


I Bp 

P+/S 


+ D, we 


V = 



1 / fiP 

a-/?\p + a 



.F 


In the latter case, 

{P + JY + (£>2 


L p^ + 2yp \ 
\ (p+y)^+w7 


• F. 


Having taken out the ‘ 1 ’ term, we need no partial fractions. 
The standard forms give us 


w = f|i — e ■>'*^cos(i)i + ^sin&»ij|. 


In the intermediate case where = iLjC, 


(p + a)' 


,.F 


/ p^ + 2ixp \ 

\ ip+ocfr 


V. 


Here the appropriate partial fractions are which 

p + a (p + a)^ 

leads to 


p 


ap 


p + a (p + a)^, 


)■ 


F 


= F{l-e-“'(l + a<)}. 

Thus, for small values of B, the condenser voltage oscillates 
about its final value, gradually settling down to it. As the 
‘damping component’ B is increased, we reach a ‘critical’ 


BO 


MOEE COMPLICATED CIECDITS 


value where the oscillations disappear. Typical examples, with 
damping half critical^ critical and twice critical, are illustrated 
in Fig. 10. 

Example 2. What alteration is made in the solution of Ex- 
ample 1, if the condenser has an initial charge to voltage 

The p-equation will now have an added Vq term, so that it will 
become 

{LCp^ + BCp+l)v= V + (LCp^ + BCp)vQ 

= (V-v,) + (LCp^+BCpi-l)v,. 

Therefore v = Vq-\ 

p^ + 

so (v — Vq) replaces the ‘v' of Example 1, and {V — v^) replaces F. 
The voltage changes from to F according to the same law by 
which it changed in Example 1 from 0 to F. 





Example 3. An impulse generator, which consists essentially of 
a condenser Ci charged to a voltage F, applies voltage suddenly 
(by the sparking over of a spark gap S) to a circuit containing 
resistances B^, B^ and the load capacitance Cg (previously dis- 
charged) (Fig. 11). Find the wave-form of voltage v across the 
load, given that 

Cl = 0-02 microfarad, = 0*001 microfarad, 

Bi = 5000 ohms, i ?2 == 1000 ohms. 

The differential equations of this circuit are most conveniently 
written in terms of v\ the voltage across at any time, as well as 
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of i;; for we know the initial values of both these quantities. If we 
wrote them in terms of only, we should need to exercise our 
minds to think out the initial value of dvjdt in terms of F. 

The current through is C^dvjdt. Passing through this 
causes a voltage drop R^C^dvjdt; therefore 

R^C^^^ + v = v. 


The current drawn from Cj_ is -C^dv'jdt. This supplies the 
current v'jR^ through and the current C^dvjdt through 

and 62* Therefore 


dv 

''‘tG 


Ri ^ dt* 


Rewriting these two equations, 


-^2^2 ^ + V = 0 , 


-ry ry dv' ry ^ dv 

B^C^— + R^C,j^+v = 0 . 


The ^>-equations are 


— v' + (R^C^p + 1) = 0, 


+ l)v'+Rfyv = RjC^V; 

for when t = 0,v' = V and v = 0. Eliminating the unwanted v', 
we find 

^ ~ {R,C,p + 1 ) {R,C,p +1)+ Rfy ■ • 

The reader will usually find it advantageous to use symbols 
until this j^-expression is reached, and then to insert numbers. 
When the numerical values are substituted, the units of (7 and 
i should be chosen so that the answer ma;Jr be in terms of a 
convenient time-unit. Thus: 

R in ohms, L in henries, G in farads, lead to t in seconds; 

R in ohms, L in millihenries, C in millifarads, lead to ^ in 
milliseconds; 

R in ohms, L in microhenries, C in microfarads, lead to t in 
microseconds. 
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In this problem, the times involved are of the order of the time of 
discharge of 0*02 microfarad through 5000 ohms — a few score 
microseconds. We therefore select the microsecond as our time 
unit. 

In numerical terms, 

'O — lOOp y 

lOOp^ -f 106p -f 1 ’ 


Y 

p2+i.o6p+o*or * 



The factors of the bottom line prove to be (p + 1*05) {p + 0-0095). 


(p 4- 1-05) (p 4- 0-0095) p 4- 1*05 "^p 4- 0-0095 
C=0, ^ = -0-96, 5 = 4-0*96. 

«; = 0 - 96 ( ^ E. \ Y 

\p 4- 0-0095 p 4- 1-05/ 

= 0-96F(e“0*oo95;_g-i.o5^^^ 

The form of this voltage is shown in Fig. 12. It rises steeply to 
its maximum valuOj then falls away more gradually to zero 
These three examples adequately iUustrate the main proposi- 
tions of the present chapter. Many more might be added but it 
seems best to defer these until after Chapter v, when a o-reat 
Simplification of th© method will be given 


then 

Hence 


CHAPTER IV 


THE RESPONSE OF A CIRCUIT TO 
A SUDDEN SHOCK 


4* 1 . Unit Function* 

A differential equation is always the expression of a process of 
change in a system. For instance, the differential equation of the 
circuit of Pig. 13, namely 


might be written 



4-^^ = F, 


dv 

It 


1 

RG 




This might be verbally expressed thus: ‘ So long as the condenser 
voltage V differs from the charging voltage F, it will increase at a 

Thus the equation focuses our atten- 
tion on the process of change. We shall 
presently find, however, that the Heavi- 
side method lends itself to a treatment in 
which the cause of the change becomes 
the centre of interest. We may introduce 
this treatment by asking ourselves how the process of change in 
the circuit of Fig. 13 is initiated; the answer is: ‘By the closing 
of the switch 8y which causes the voltage across the terminals AB 
to rise suddenly from zero to F.’ This sudden change of voltage 
gives the circuit a shock, to which its response is expressed by 
the differential equation. 

If F had been, for example, an alternating voltage, the shook 
to the circuit would have been different, and of course its response 
also. However, we select one form of shock as our ‘standard 
shock’, and build up our theory from it; just as alternating 
current theory is built upon the sine wave, and waves of more 
irregular shape are treated as a combination of sine components. 
The standard shock is a sudden rise (of voltage, current, etc.) 
from zero to 1, taking place (unless otherwise stated) at time 
sg = 0. Before, the quantity is supposed to have the constant 

XET 


B 




Fig. 13. 
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value 0; afterwards, the constant value 1 (Fig. 14). This standard 
shock is called unit function, and is commonly written 1, or 
sometimes S(i5) : H suggesting its originator, Heaviside. The latter 
notation has the advantage that it lends 

itself to the expression of a unit function 

whose leap occurs at time tQ, for this may 1 

he written H{t — 1^ \ but the notation 1 is — ^ 

more compact for ordinary purposes, if 

this refinement is not required. 

A leap in voltage from 0 to F, such as occurs across terminals 
AB of Fig. 13, will therefore be written F x 1 , or F . 1 . The closing 
of a switch in a d.c. circuit is in fact one of the commonest ways 
of applying unit function voltage to a circuit, while the opening of 
such a switch frequently causes a current impulse of unit function 
form. These points will be clarified later — ^they are mentioned 
here to show the fundamental importance of the unit function. 

4*2. Operations upon Unit Function. 

Suppose that a circuit is shocked by unit function voltage F. 1, 
and that the voltage v in some portion of the circuit is obtained 
in the form 

^ = F(p).(F.l). 

We already know how to interpret F{'p) operating on F ; what if 
F{f) operates on F. 1 ? 

The answer to this question takes us right back to our defini- 
tion of p, which was 



It will be seen that p only takes account of times subsequent to 
i = 0; therefore the form of the solution is just the same, whether 
F(p) operates on 1 or on 1 . These two statements crystallize the 
two ways of looking at the matter. 

(1) ' F(p) operating on 1 leads to some function of say W{t ) ; 
but this must only be used for positive values of negative values 
are of no significance.’ 

(2) ‘ F{p) operating on 1 leads to W{t) x 1, which is zero when 
t is less than 0, and equal to W{t) when t is greater than 0.’ 

For example, the circuit of Fig. 13 leads to the solution 

t? = F(l-~-e“*«^)xl. (Compare §2*3, with ^7o = 0.) 
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TMs is illustrated in the full curve of Fig. 15. The solution in § 2*3 
would have led to the dotted curve, if we had applied it to the era 
before the closing of the switch; but that would have been 
foolish, because it is by closing the switch that the whole process 
is started. The two methods are exactly equivalent if the earher is 
seasoned with common sense. Consequently, in the equation 

F(p).t^W{t)xl, 

the 1 on the right-hand side is nearly always omitted. 

Where attention is directed to the cause of change, not the 
process, the reader will find that operations on 1 seem more 
natural than operations on 1 ; this, however, is just a philosophical 
distinction, which does not affect the final result. 



4-3. Example. Response of a Circuit to Unit Function 
Voltage. 

A transformer is connected to an overhead line through a cable 
of capacitance 0*01 microfarad, and a reactor of inductance 40 
millihenries; the transformer itself can be regarded as an in- 
ductance of 5 millihenries. A hghtning stroke causes on the line 
a very rapid voltage rise of magnitude F (Fig. 16). Find what 
happens at the transformer terminal. 

Let the voltage across the transformer winding be v. 

Then the current in C = Cdvjdt, 

Therefore the rate of change of current in (7 = 0 d^ldf. 

Also the rate of change of current in L ^ vjL. 

Therefore the rate of change of current in I = the sum of these, 

= Cdhldt^-^v/L. 

Therefore the voltage across ? = IGd^^vIdt^ + lv/Ly 
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and the differential equation is 



with the initial conditions v = 0, dvjdt = 0. The ^-equation is 


Substituting iy = 5xl0® microhenries, ^ = 40x10^ micro 
henries, C = 10“^ microfarad, 


whence 


the multiplying 1 being omitted. This voltage wave is illustrated 
in Fig. 17. 

Though simple, this example 
containstwofeaturesofinterest. 

(1) The unit function voltage v/ / \ / 

wave is employed as an aunroxi- ® ' / \ / 


Microseconds 
Fig. 17, 
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supposed source of unit function voltage, wMcli must have an 
effective impedance of zero, because the voltage which it pro- 
duces is unaffected by the load. 

4*4. The Principle of Superposition. 

We tend to think of the voltage in a circuit as causing the current, 
but there is no logical reason why we should; often it is useful to 
think of the current as causing the voltage. It may even be more 
convenient to look upon some other quantity — such as the charge 
on a condenser, or the magnetic flux in an inductance — as the 
prime cause of the circuit’s behaviour. Let us embrace all these 
possibilities with the non-committal title ‘Stimulus’, and the 
symbol U, V will generally be a function of time; thus, the applied 
stimulus may be a voltage of any wave-form, applied between two 
chosen points in the circuit. The term ‘shock’, already used, 
means ‘a stimulus suddenly applied’; the term ‘stimulus’ in- 
cludes ‘shocks’ also. 

Suppose the circuit is initially dead; that is, no condenser is 
charged, no inductance is carrying current. Let a stimulus be 
applied to the circuit. In each branch, voltage, current, charge 
or flux will now come into being; but we are probabty interested 
in one property of one branch, so we will fix attention upon it, and 
call it the ‘ Response ’ u. Thus, we may be chiefly interested in the 
voltage across a particular condenser. 

The Principle of Superposition may be stated thus: If the 
stimulus U applied to a circuit causes a response u in a particular 
branch, and if the stimulus U' causes in the same branch a re- 
sponse u', then U and U' acting together will produce a response 
u + u'. (Note that there is no need for U and V to act at the same 
points in the circuit.) This principle is an immediate consequence 
of our basic assumption (§1*1) that each circuit element is linear. 
For example, 

Current i, in an inductance L, produces a voltage 
V = Ldijdt; 

Current in the same inductance, produces a voltage 
v' = Ldi'jdt; 

so {i+i') produces a voltage L{dildt+di'jdt), oi {v + v'). This 
would not be true, for example, in a saturated iron-cored in- 
ductance. 
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4-5. Consequences of tlie Principle of Superposition. 

The following consequences scarcely need justification : 

(1) If a stimulus is applied to a circuit which is already alive, 
without the characteristics of the circuit being otherwise altered, 
the resulting voltages, currents, etc. may be calculated by adding 
together 

(а) the voltages, etc., which would have been there if the 
stimulus had never been applied, and 

(б) the voltages, etc., which the stimulus would have caused 
if the circuit had been dead. 

The trouble about this is that most ways of applying a stimulus 
to a circuit (like closing a switch) do alter the circuit character- 
istics at the same time, and this alteration has to be allowed for. 
We shall see how to do this in the next chapter. An example of a 
stimulus which may be applied with only slight alteration in the 
circuit characteristics is a voltage applied by exciting an alter- 
nator, already connected in circuit. 

For the rest of this chapter, we shall assume that the circuit is 
initially dead. The results of any initial 'aliveness’ can be dealt 
with separately. 

(2) If the stimulus applied to a circuit consists of a succession 
of shocks, the response to the stimulus may be obtained by adding 
together the responses to the shocks. 

4*6 . Shocks which are not of Unit Function Shape . DuhameFs 
Integral. 

This last statement shows us how to determine the response of a 
circuit to a shock of any shape, applied at time t = 0, w^hen once 
we know its response to a unit function shock; for the general 
shock may be treated as a succession of shocks each of unit 
function form (Fig. 18). 

Suppose we know that unit function shock evokes a response 
^ 1 , which is a function of t and may therefore be written Ui(t), 
What happens when we apply a shock of form U (t), at time t = 0? 

Right at the start, the circuit receives a shock 17(0). 1, which 
evokes a response 

U(0).Uj(t), 

After that we get the eiOfect of a succession of little jumps of unit 
function shape. If we desire to know what the circuit is doing at 
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time t, we must take account of all the little jumps which are 
engendered up to that time. Let us focus attention on the one 
which starts at time r, and let us suppose that a time Sr elapses 
before the next one. Then the height of that particular jump must 
be U'{r)xSr, where U'{t) denotes the value of dU/dt — the 
gradient of the Z7-curve — at that particular instant. 



Thus at time r a shock U'{t) Sr,li& applied. When time t comes 
round, this shock has been eifective for a time (t—r), and so has 
caused a response 

U'{r)ui{t — T)ST, 

Consequently the response of the circuit to the whole succession 
of jum]3S up to time t is 

2 U'{r)u^{t-r)ST, 

T = 0 

and proceeding to the limit where St becomes indefinitely small, 
we find that the shock U{t) causes a response u{t), given by 

u{t) = U{0)u^{t)+ j D'{T)u^{t-T)dT. 

This is called Duhamel’s integral. Notice that, in working out 
the integral, t is the variable; t is treated as a constant. 

By integrating by parts and by other elementary means, we 
may prove that Duhamel’s integral can be written in the following 
alternative ways: 

(1) U{t) = U{0)u^it)+^ U'{T)Ui{t-T)dT. 

(2) u{t) = u-^{0) U (t) + j u[{T)U{t- T)dT. 

(3) «(i) = 17{0)%(i)+J^%(T)t7'(«-T)dT. 
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(4) 

rt 

u{t) = %(0) ?7(^) -j- U{t) u[\ 
Jo 

(5) 

II 

|* f7(T) — t) c^rj 

(6) 

II 

1* %(t) V{t — r)dT^ 


Comparing (1) with (2), (3) with (4) or (5) with (6), we see that 
U and Ui are interchangeable. But the presence of U{Q) in the 
first term is always followed by U' in the integral, while Uj^(0) in 
the first term is followed by in the integral; to notice this may 
assist the reader’s memory. 

4-7. Examples in the Use of DnhameFs Integral. 

Example 1 . In the circuit of Fig. 1 6, the lightning stroke, instead 
of causing a voltage change of unit function form, causes a voltage 
wave 

F(e”^ — (Fig. 19), 

where a = b — f . Find the response voltage across the trans- 
former. 

The solution already obtained 
shows that the response to unit 
function voltage is given by 


%(^) - — 


o-zv 


(o = Also 

,-U\ 



where k = 

Uit) = We ob- 

serve that will be simplified 
by differentiation, while U lends 
itself more kindly than to 
having t replaced by (i — r). (% 
would split into a sine and a cosine term, whereas V remains single.) 
We therefore select the second form of DuhameFs integral, and 
obtain 

u{t) = 0 + J ko) sin WT . dr 

= sin sin cot)} dr, 

1 

Now e«^sm<or(^r-~ 2 ^{cxe«^sinco^ + 

*/ 0 w 1 CO 
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Therefore 
u(t) = Vkci) 

= Vkoj 


a sin (j)t — (i) cos (jjt + (j)e~^^ b sin (ot — o) cos 4- m~ 


b — a 




- sin — <9) + - 


(jt) 


(t) 


,V[(«^ + io^) {b^ + OJ^)] ' a^ + co^ ^ 
where tan d — (j)(a + 6)/(a6 — 

With numerical values inserted, this becomes 


52_|_^2 


o-U 


Mt) = sin 0-865) + ^e- 




1681 '' 


« -21/31. 


The form of this voltage is shown in Fig. 20. 



Example 2. DnhameFs integral deals very neatly with wave 
shapes which, though capable of representation by mathematical 
functions, have points of discontinuity where the function must be 
changed. Examples of such waves are the saw-toothed voltage 
forms employed in cathode-ray tube sweep circuits (Fig. 21), and 



the voltage forms produced by rectifiers (Fig. 22). The following 
is an instance of the latter type* 

The simple smoothing circuit (r, C) (Fig. 23) is supplied from 
a mechanical rectifier which generates half-waves of form F sin Qt, 
and is loaded with a resistance J2. Find (in the steady state) the 
wave -form of voltage across the load, 

(The mechanical rectifier is specified so as to avoid the added 
difficulty of dealing with a rectifier hke a thermionic valve, which 
cannot pass a reverse current.) 
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If V is the voltage across the condenser, the current in C is 
Cdv/dt, and the current in M is vjB, Therefore the voltage 
across r is 


i^dv v\ 


and if unit function voltage is applied to the circuit, the dif- 
ferential equation is 


„dv 


The ^-equation is rC{p + cc)v 
to this is 


B -j- T 

1 , where a = - 77 ^ . The solution 


CBt' 


V = 


1 - 


and this is the function which we called 

We now deduce from this the effect of shocking the circuit with 
a single half-sine-wave Fsini3^ (Fig. 24), which, since it only 
lasts from ^ = 0 to 7r/i2, should be written 


r "iTr/n 

U{t) = I Fsini3^ 


We will calculate the response for times greater than njQ — that 
is, for times after the half-sine-wave has come and gone. 

We note that u{{t) is simpler than and also that u[{t) lends 

itself better than U{t) to having t replaced by {t — t). We therefore 
select the fourth form of DuhameFs integral, and obtain 


^(or u(t)) ■ 


J‘[Fsmfc] 

i ^71 

Jo 


0+ 

Tc 


IT IQ p—CX,{i—T) 


rC 


■dr 


BmQrdr, 


since the period from time ttJQ to t contributes nothing to the 

integral, ^ 

[a 


rC ljx^+£2^ 


^ 0,7 "Itt/D 

(a sin — i2 cos J0 t) 

Jo 


-cd 


rG{a?+Q^) 

If we measure t from the end of the half-cycle, we get 

Fi3(l - 1 - 6-“’^/^) 


rC{o!?-{-Q^) 


-cd 
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The contribution of the previous half-cycle is got by increasing 
t by ttIQ-, so it is 

^ i Q). 

rO(cc^ + Q^) ® 

the one before that contributes 
Fi3(l + 


V = 


TC{a^-\-Q^) 




and so on. Therefore all half-cycles before the present one 
contribute 

V Ql \ 4 - p-^^f^\ 

^ '^-ai ^ 1 p—cc7TjD ^ p^-2a7rl£) J 


rC{(z^ + Q^) 
VD 1 + 


rC{(X^ + Q^) 


)—oU 


provided t is taken as being less than irjQ. 



Fig. 24. 



Finally, the present half-cycle contributes 

Ye-ccin 

V = — -pr- «“’■ sinl^rdr 

fC Jo 

Yp—oii r paT 

F 

= ^ ^2) Qt+Qe-’^). 

Therefore the total effect of all the half-cycles is 

Fi3 (O^ • ru 2 A 

^ ~ r C(a2 + Q^) \i3 - cos -F ^ e “ j . 

The form of this voltage is shown in Fig. 25. This method 'was 
suggested to the writer by his colleague, Mr D. J. Mynall. 
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The alternative method of solving this problem would be as 
follows; 

(а) Assume a value for the voltage at the beginning of the 
cycle. 

(б) From the differential equation of the circuit, write down 
the p-equation utilizing this initial value v^. 

(c) Solve the p-equation, and from the solution find the value 
of V at time njQ. 

(d) Equate the initial and final values of v, thereby ■marlring 
the fact that the circuit is in a steady condition. This equation 
will give the unknown Vq. 


4-8. Note on the Application of Duhamel’s Integral to the 
Aftermath of a Pulse. 

A shock which terminates after a finite time, like the half-sine- 
wave of Fig. 24, may conveniently 
be termed a ‘pulse’. This particular i 
pulse is made simpler to deal with by 
having no ‘step’ either at beginning 
or end; if it had terminated with a 
step, like the pulse shown in Fig. 26, 
this would have to be allowed for, 
thus: 



A 

Fig. 26. 


The initial step contributes 17(0) %(f) to the integral. 

The intermediate part contributes J*' TJ'{r) u^{t - r) dr. 

The final step contributes - 

The sum of these three may readily be transformed into the single 
expression , 

U{T)u[{t — T)dT, 

which resenibles the fourth form of Duhamel’s integral, except 
that the initial term is missing. This is the best form of the 
integral to use, when calculating the after-effects of a pulse upon 

,a'" circuit. ' 
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4-9. Gomparison of Methods for dealing with Shocks which 
are not of Unit Function Shape. 

The reader will probably have recognized that, when the dif- 
ferential equation of a circuit is written down in some such 
form as 


a 


d^u 


+ b 


du 

dt 


+ CU =f{t), 


the f(t) on the right-hand side is closely related to the stimulus 
which is applied to the circuit. (See, for instance, Example 2 in 
§ 2*8.) Therefore, contrary to the canons of economical teaching, 
we have learned three methods of dealing with shocks which are 
not of unit function shape: 

(a) In Chapters ii and iii it was suggested that operational 
expressions like 


1 

p + a 


m 


should be obtained, and interpreted by the rule 

dt. 

(6) In Chapter m f{t) was to be written in the form <p{p),li 
and combined with the other p-expressions, the whole lot then 
being split up into partial fractions and interpreted. 

(c) In Chapter iv we have proposed solving the problem for 
a unit function shock, and then deducing the effect of the more 
complex shape by the use of Duhamel’s integral. 

I suggest that the choice between these should be made in the 
following way: 

(a) is really equivalent to (c), and may be considered to have 
been superseded by it, except perhaps for simple problems; 
(5) involves algebraic processes, which are more elementary and 
less liable to error than the integration involved in (c). Therefore, 

When the stimulus is a simple continuous function, use (5); 

When the stimulus is capable of mathematical representation 
but has discontinuities, or when no mathematical form is known 
for it, use (c). (If the stimulus is given as a curve of no known 
mathematical form, Duhamers integral can be used in conjunc- 
tion with some method of numerical integration; see § 7’4.) 
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440. Tlie ^Hammer-Blow’ Shock. 

Besides unit function, there is another standard shock which is 
worth mentioning: a pulse of extremely high value and short 
duration, like the hlo'w of a hammer. Such shocks are well known 
in dynamics, under the title of impulses 5 unfortunately this 
word is used rather more vaguely in electrical circuit work, so it 
has seemed best to avoid it here. A dynamical impulse is an 
infinite force lasting for an infinitesimal time, but having a time- 
integral which is finite. The corresponding electrical impulse will 
therefore be an infinite current lasting for an infinitesimal time, 
but transporting a finite charge; or an infinite voltage lasting for 





t-i- 
Fig. 27. 

an infinitesimal time, but associated with a finite change of 
magnetic flux. Infinite voltages and currents do not occur in 
nature, yet they may form a convenient approximation to the 
truth; for if the duration of a stimulus is small compared with the 
natural periods of a circuit, there is little error in treating it as a 
'hammer-blow’. 

Fig. 27 indicates how to find the operational form of such a 
shock. The pulses U{t) all enclose an area 1; the integral curves 
shown in the right-hand column all rise to a height 1 . As the 
pulses get shorter, the fronts of the integral curves get steeper; 
plainly, when the height becomes infinite, 


-U{t)dt 
% 0 


1, 


or 
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-* 7 ( 0 = 1 , 

P 

or U{t)=p.l. 

This is the operational form of the standard ‘hammer-blow’ 
shock. 

411. Example. ‘Hammer-Blow’ Shock. 

The circuit of Fig. 16 is struck by a lightning stroke, which causes 
a steep voltage rise F ; but after 1 microsecond this is brought 
equally abruptly back to zero by an insulator flashover. Find the 
voltage at the transformer terminal. 

The natural period of the circuit is 4077-/3, or about 42, micro- 
seconds, which is long compared with the duration of the pulse; 
so we treat the pulse as a ‘hammer-blow’ of the same area, 
namely V volt-microseconds. 

The differential equation is now 

I rr ^ 

^0^, + j^v + v=Vp.l, 

so the j3-equation is 



Substituting the values of L, I and C, we get 

= p JL 

^ F^+(^)‘400' 

V , 


Comparing this with the expression obtained for unit function 
shock, we see that the maximum amplitude of v is here only 
F/60 instead of 2V {9 . Of course the 'hammer-blow ' shock cannot 
give just the same wave-form as the 'chopped unit function’ for 
the first microsecond; the latter must start in the mode shown in 
Fig. 17, for how is the transformer to know that the voltage is 
going to be chopped off, until the chopping happens? But the 
expression just found is a very good approximation to the subse- 
quent oscillation of 
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IMPEDANCE AND ADMITTANCE OPERATORS 


5T. The Analogy between Alternating and Transient Con- 
ditions. 

When an alternating voltage V is applied to a circuit, whose 
admittance is known to be 

G+jB, 

the components of the current taken by the circuit are given by 

{G+jB)xV. 

Furthermore, if the circuit has several branches, the current in a 
branch n will be proportional to F, and may therefore be written 

(G^+jB,)xV, 

where (Gn+j^n) is another admittance associated with the cir- 
cuit, sometimes called the ‘transfer admittance ’ of the branch n. 
Thus, in Fig. 28, the admittance of the whole circuit is 

Gjo) 1 

rCjo) -i- 1 Lj(t) + B ’ 

while the transfer admittance of the condenser branch is 


Cj(t) 

rOjo) + 1 ■ 


Now suppose that V, instead of being an alternating voltage, 

is a voltage shock applied to the 

terminals of the circuit, which was ^ i 
previously dead. From the pre- [frequency ctf/2?r] 
ceding chapters, we know that the | 

current which flows into the circuit Fig. 

can be obtained by operatingupon V with somep-expression, while 
the current in any particular branch will be obtained by operating 
upon F with some other p-expression. In this circuit, for instance, 



Current 


and 


(rCp -hi + i?) • 


[rCp - 1-1 ip - 
Current in condenser = ■ 1 ). 


FhasbeenwrittenF(i),toremindusthatitmayhaveanywave-form. 
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There is evidently a close analogy between the a.c. and the 
transient expressions. The currents are obtained from the voltage 
function V{t) . 1, by operating upon it with p-expressions known 
as admittance operators. In exactly the same way, the voltages 
associated with an apphed current shock are obtained by 
operating upon it with impedance operators, corresponding to 
the A.c. equation 

F = (i2+jZ)x/. 

kloreover, the voltages associated with a voltage shock are ob- 
tained by using operators which are reaUy the ratio of two 
impedance operators, just as in Fig. 28, 

Voltage across L = ^ 

Lj(j) + R 

This little example will have suggested to the reader a fact full 
of significance: namely, that the impedance and admittance 
operators have just the same form as the a.c. impedances and ad- 
mittances, except that jo is replaced by p. That this must be so 
under certain conditions, we shall now proceed to prove. 


5-2. Impedance and Admittance Operators in Circuits 
Initially Dead. 

We shall first consider circuits which are dead, up to the instant 
of apphcation of the shock to which the response is desired; by 
this, we mean that no inductance is carrying current, no con- 
denser is charged. Another way of saying this would be to state 
that the circuit has, initially, no stored energy. 

At t = 0 the shock begins. On each condenser C, the ensuing 
voltage and current are connected by the relationship 


if the positive direction of current is taken to be that associated 
with an increase in voltage. Operate on both sides of this equation 
with the operator 


We get 


1 

V 

1 . 


= Cv, 


CET 


4 
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since = 0 when = 0. In other words. 


On each inductance L, the voltage and current are connected by 


if the positive direction of voltage is taken to be that associated 
with an increase in current. Therefore 


since i = 0 when ^ = 0. In other words (algebraically) 


Exactly similar reasoning applies to each mutual inductance M 
here we shall get 


where 1 and 2 of course denote 'primary’ and 'secondary’. 

We may therefore formulate the following rules for finding the 
response when a dead circuit is activated by a shock. 

(1) Write down the relationship between the required re- 
sponse u and the known shock C7(^) . 1, treating the circuit Mke 
an A.c. circuit in which each condenser has impedance IjOp and 
each inductance has impedance Lp or Mp. 

(2) Interpret the resulting expression, whose form will be 

u = F(:p).{U(t).i}, 

by carrying out the operations denoted hj F{p). 


5-3. Examples. Impedance Operators in Circuits Initially 
Dead. 

Fxample 1. The circuit of Fig. 29 is dead until the switch /S is 
closed, connecting it to a d.c. supply of voltage F. Find the 
voltage V across the condenser. 
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The rule jiist given leads to the result 


or 


V __ l/Cp 

V ^ L^+B + (iiGiy 

LCp^ + RGp + V^^-^'^' 


The rest of the solution is just the same as is given in §3T1, 
Example 1 ; but see how much more quickly the p-equation has 
been obtained ! 



Fig. 29. 



Fig. 30. 


Example 2. In a certain oscillograph circuit it is necessary to 
step down a voltage of half-cosine-wave form (Fig. 30) by a con- 
denser potential divider G. When an inductance L is connected 
in parallel with C (Fig. 31), it is found that the stepped-down 
voltage is distorted in the manner shown in Fig. 32. The second 
positive peak is undesirable. Account for this peak, and suggest 
a cure. 


“ Vco/j 

L.. 


Fig. 31. Fig. 32* 

We take the input voltage as {VoosQt)A, up to the time 
t = ttIQ. The output voltage will be v, as usual. 

The operational impedance of the circuit of C in parallel with 
liis 

1 Lp 

Op + {ljLpy ~ LCp^+ V 
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.{(F COS 1} 


Lp 

Therefore v = — . {( F cos . 1} 

LO 

L{G + Co) + V 

We shall interpret this by writing (cosfi^).l in its operational 
foriEj ^2 _j_ jQ 2 * ^ * ^^CrCOVCIj let ^ ^ ^ (0l2i7T is ijlio 

natural frequency of the circuit. Then 


C+Co'ip^+a)^){p^+Q^) 

- "^^0 / j 5 !£!_ "V \ ^ 

(C + Go){Q^-o)^)\p^ + Q^ p^ + ojy-^ 

_ VCo /i3^coS;Qi — o^coswA 

~ c+^\ )■ 

During the interval from i = 0 to TrjQ, cos Qt diminishes from 
+ 1 to - 1. If w is greater than Q, eosot will diminish from + 1 
to - 1 and then increase again, so that it is quite possible for the 
whole bracketed expression to be positive in the neighbourhood 
of t = nlQ. For example, Fig. 32 has been drawn with co = BQj2. 
On the other hand, if w is less than i3, the bracketed expression 
cannot become positive again; therefore a remedy for the second 
positive peak is to make a> less than D— that is, the natural 
frequency of the circuit should be lower than the frequency of 
the applied voltage. 

5‘4. Shock due to Closing or Opening a Switch. 

Probably most transient conditions in engineering circuits are 
brought about by closing or opening a switch, or by some equi- 
valent operation. Closing or opening a switch has two effects 
upon a circuit: 

(1) It applies a voltage or current shock to the circuit, 

(2) - By altering the interconnection of the circuit elements, it 
changes the operational impedance of the circuit. 

Consider first the effect of closing a switch. While it is stiff open, 
there will be in general a voltage across its terminals which may 
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vary with time in any manner; we will call this voltage /(i) 
(Fig. 33). In the various arms of the circuit, there will be voltages 
Vi and currents associated with the ‘open’ condition of the 
switch; these may also vary with time in any manner. 



Now suppose we have another exactly similar circuit, initially 
dead as regaids both inductances, condensers and generators 
but having instead of the switch a generator of zero impedance" 
which generates a voltage -/(«).! (Fig. 34). This generator 
applies to the circuit a voltage shock which leads to response 
voltages and currents 




Superimpose these two systems of voltage and current. Across 
the switch we have voltage f{t) before t = 0, and zero voltage 
afterwards — just as if it had been closed. Closing the switch is 
therefore equivalent to applying a voltage shock -/(«) . 1 ; and 
the voltages and currents which result from closing the switch 
are + ii+ig. In other words: 

The changes produced in a circuit by closing a switch are the 
same as would be produced by injecting across the switch 
terminals into a dead circuit a voltage equal and opposite to that 
across the open switch. 

In exactly the same way, we may prove that 

The changes produced in a circuit by opening a switch are the 
same as would be produced by injecting across the switch termi- 
nals into a dead circuit a current equal and opposite to that 
through the closed switch. 
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It should be added that, in calculating the response of the 
circuit to a voltage or current shock —f(t). 1, all voltage gene- 
rators in the circuit must be treated as being of zero impedance, 
all current generators of infinite impedance. 

5*5. Impedance and Admittance Operators in Circuits not 
Initially Dead. 

The whole basis of impedance operators, as so far described, rests 
upon the assumption that the circuit is initially dead. This is 
obviously a grave limitation upon the method. Of course the 
methods of Chapters ii and m gave us a way of treating circuits 
not initially dead, for in such a p -equation as 

{ap^ ‘j-bp-{-c)u — f{t) -f {ap^ -i- bp) Uq H- ap) 

the Uq and terms represent the initial state. But that involves 
writing down the differential equation, and jettisoning the pretty 
simplicity of the impedance operator method. Can we not retain 
this simplicity, even when the circuit is initially alive? 

Broadly speaking, this end can be attained in two ways. One 
way is to split up the problem into parts which each involve a 
shock on a dead circuit; the principle of superposition justifies 
this. The other is a rule of thumb method for writing down the 
; 2 ?-equation for any initial conditions. 

5*5 (a). Superposition Method. 

In § 5*4 we saw that the change produced in a circuit by closing 
(or opening) a switch could be calculated by treating the circuit 
as dead. Adding this change to the conditions which would have 
existed if the switch had not been closed (or opened), we obtain 
the actual state which results from the switching operation. 

Now in many cases the conditions, which would have existed 
if the switch had not been operated, are extremely simple. Often, 
for instance, those conditions are steady — perhaps there may be 
a condenser which remains fully charged until a switch is closed. 
In such a case the only real problem is to calculate the change due 
to the switching; and this may be done by direct application of 
impedance operators. 

Example I, In the impulse generator circuit of Fig, 35 the con- 
denser Oj charged to voltage F is suddenly connected to the rest 
of the circuit by the sparking over of the spark gap 8, Find the 
voltage V across the load condenser C 2 . 
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Before tlie sparkover of the gap, terminal A is at potential F 
with, reference to terminal B. The sparkover applies a shock 
— F. 1 to A relatively to B, or h-F. 1 to B relatively to A. 
So the changes produced in the circuit of Fig. 35, when S becomes 
conducting, are the same as the response of the circuit of Fig. 36 
(initially dead) to a shock F . 1 . 


Vq 


T 


i-p 

A B 




'/c, 


■ 2 /> 


Fig. 35. 


B 



Fig. 36. 


The operational impedance of Cg with R in parallel is 


1 B 

C'2p + (l/i?)’ ~ RG^ + V 

R 

Therefore v — — . F. 1 

^ Y 1 

{LC^'^+1){BG^P^+1) + RCyp ' 


If 8 had never been closed, v would have been persistently zero; 
therefore there is nothing to add to this expression for v. The 
factors of the bottom line in this p-expression can take different 
forms according to the relative values of i?, and Cg, so this 
is the point at which numbers should be inserted. 

To illustrate the necessity, sometimes, for adding the condi- 
tions which would have existed if the switch had not been closed, 
let us calculate the voltage on (7^. From Fig. 36, the change in this 
voltage is 

1 




B r 1 


V.l 


^ y 1 

~ {LC^^+l)iRO^+l) + RGyp 
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The polarity of this change is such as to make A the 'negative^ 
terminal j which is opposite to the initial polarity. Therefore ^ 
superimposing the change upon the initial conditions, we find 


Voltage on — 


+ 1 y I 

(LG^^+l)(BC.,p + l) + RC^p ' 


Example 2. A generator field of resistance R and inductance L 
is drawing a steady current from an.o. voltage source F, when the 
contact ^ of a Tirrill regulator is opened, inserting into the 
circuit the resistance r and condenser G (Fig. 37). Find the 
voltage V which appears across the regulator contacts. 



Initially, the contacts were passing a steady current F/i?. The 
effect of opening them is to apply a current shock (F/i2) . 1 across 
them; and, as has been explained already, the voltage source 
must he treated as having zero impedance. Fig. 38 therefore 
shows the equivalent problem. 

The three current paths are in parallel. Their operational 
impedance is 

1 ^ riLp + R) 

^ I ^ ~’rLCp^+{ErC + L)p + {R + ry 

^ r Lp + R 

Therefore the voltage across the contacts is 

r{Lp + R) V 

rLGp^ + {RrC + L)p + (R + ryR 

P + {ii/L) X 1 
{p-\-a)ip+j3yRC ' 

where a+/3 = {RrC + L)lrLG, aji = (R + r)lrLG. 

(This would be the convenient form to use if the factors of the 
bottom line were real.) 
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Taking partial fractions 


we find 


p + {E/L) _ Ap Bp 
{p + cc){p + ^) p-i^a 


Therefort 


{BIL)-a 
a{a—/3) ’ 


{RIL)-I3 R RfC 

■ ’ Loc^~ R + f- 


RG\a{a-p)^ /?(a-y?) 




V R + r^LC(a-jff)[{cc R, 






There is nothing to add to this expression, for v would have re- 
mained zero if the contacts had not been opened. 


5*5 (6). Rule of Thumb Method* 

Where the state of affairs before the switching is less simple, and 
in all cases of doubt, it is useful to be able to fail back upon a rule 
of thumb method for obtaining the ^-equation. 

Suppose a cffcuit has n elements, that is, the total number of 
resistances, inductances and condensers is n, each mutual in- 
ductance counting as two. We know all about its behaviour, if we 
know the voltage v and current fin each element. To find these 
2n unknowns, we need 2n differential equations, n of them are 
provided by the relation between voltage and current in the 
individual elements: v = Bi, v = Ldijdt, etc. The other n are 
provided by Kirchhoff’s well-known laws of circuit behaviour, 
namely 

(1) The sum of the voltage drops round any closed circuit is 
equal to the sum of the e.m.f.’s in that circuit. 

(2) The sum of the currents arriving at any point is zero. 

To write down the ^-equations of a circuit under any initial con- 
ditions, all we need is to know the jp-relation between the voltage 
and current in each element. 

The initial state of the circuit is completely specified if we know 
the current in every inductance, the primary and secondary 
currents in every mutual inductance, and the voltage on every 
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condenser. Suppose we know that an inductance L has initially a 
current The voltage across it is given by 


Therefore 

or (algebraically). 

This is the ^-relation between voltage and current in an in- 
ductance with initial current Similarly, in a mutual inductance 

= Mp(ii - ijo) + L^pii^ - iao). 

And in a condenser, 

i = Cpiv — Vo). 

We can therefore formulate the following rules for writing 
down the ^-equation under any initial conditions : 

(1) Take as separate variables the current in each inductance, 
the two currents in each mutual inductance, and the voltage on 
each condenser. 

(2) Assuming the above relations between voltage and current 
in the circuit elements, write down the Kirchhoff equations for 
voltage and current balance. (The switching operation, or what- 
ever initiates the disturbance, must be regarded as having taken 
place if the circuit connections are affected thereby.) 

(3) Eliminate the unwanted variables, thereby obtaining a 
jp-expression for the one which is required. 

This method directs attention to the process rather than to 
the cause of change, so, as suggested in §4*2, the symbol 1 does 
not appear. 

Example 1. Let us apply this method to the problem of Fig. 37. 
The variables shall be the current in L (initial value VIE), and the 
voltage across <7 (initial value 0). 

After the contacts S have been opened, we have the voltage 
equation 


^ di 

Jr 

V = Lp{i-io). 
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since the voltage on L is Lp(i—ig) ; or 

v + {Lp + R)i = V+^V; 

It 


and also tlie current equation 


Gpv + - = i, 
r 




Eliminating i between these equations, we find 


(L^) + i?)^C'i9+ij + ljv= (- 


+ 1 F, 


>nt. 1 


, r{Lp + R) V 

whence v YLCp^-^{RrC^L)p^{R^r)' R’ 

just as before. 

Example 2. In a certain oscillograph circuit, an operation 

equivalent to closing a switch 8 

connects the circuit shown in ^ ^ . 

Fig. 39 to an alternating voltage | ^ <L 1 

source at the instant of voltage ^ a > ^2 

maximum. When 8 is closed, a ^ ^ ‘ ' 

current ^20 is still flowing in the 

secondary circuit, the residue of a previous cycle. Find the 
changes in secondary current afterwards. 

By the rules given, the primary voltage of the mutual induc- 
tance is 

the secondary voltage, 

Mpi^ + L^pii^-i^o). 


Therefore 


{Ljp-{- Bi)ii-\-Mp{;i2—i^o) = Fcosi3^ 
Mpi-^ + L^p{i^-i^Q)^R>^i^ = 0 ; 

(ijj) + i?i) 4* Mpi^ = V cos Qt 4- Mpi^^, 
Mpi^ 4- {L^p + i?2) H = 


Therefore 


{(LiP 4- i?x) (4 p + ^ 2 ) - ^2 

= — MpV cos Qt 4 - {(^ 1^2 + -^1-^21^} ^ 20 * 
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This may be written 


{LiLi - if 2) + a) (p + /?) 

^ ~ ^ - -3^^) (P® + 7J)) ^20= 

where 


oc+P = 


B^L^+R^L^ 

L^L^-M^ 


5 £X^ 


-^1-^2 

ilia -if 2’ 


7 = 


Bih 


Hence 

« T7 I i>^ + 7i> • 

^ Llig — Jf2(j) + a) (^+^)(^2^^2)- (i? + a) (^5 +/?) ’ 


Each term can be interpreted by splitting into partial fractions, 
as usual. 

The reader may find it instructive to attempt this problem by 
the superposition method. The voltage, which would have existed 
across the switch if it had not been closed, is complicated by a 
voltage induced by the dying current whose instantaneous value 
was iao- The method just given takes account of all this without 
effort. 


CHAPTER VI 


RESONANCE, DAMPING AND STABILITY 

6-1. Deductions from the />-Equation. 

In the preceding chapters we have bent our energies to finding 
a complete expression for the response of the circuit as a function 
of time. But in many problems we may be content with less t,ba.n 
this; we may wish to know merely the kind of response which the 
circuit will make, or what changes must be made in the circuit to 
bring about a response of a more desirable kind. Such information 
as this can often be deduced directly from the ^ 9 -equation, without 
going through the ensuing process of interpretation. 

When we obtain a p-equation such as 


u = 


yp 


. 1 , 


(p+(x){p+^y 

the type of solution is controlled by the factors of the bottom line 
of the j5-expression; for these determine the partial fractions into 
which this j9-expression shall be split up. In the case just quoted, 
these are 

p p 


u 


j 

ct-y\p+^ P + Ol.)' ’ 


leading to 


u 


y 


a—/? 




Differences in the top line lead to different solutions of the same 
tjrpe; thus 


u = 


yp‘‘ 


gives 


{p + oc){p+^) 

y < 


.1 




a different combination of the same exponentials. 
If the j3-eqiiation were 


it would be necessary to express /(^) in the form ^(p). 1, before 
examining the p-expressions. That part of the solution which 
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depends upon and will have its type controlled by the 
factors of {ap‘^ + hp-^c). The part depending upon /(^ will in 
general contain these terms too, but it will also contain other 
terms corresponding to the factors of the bottom line of 

To determine the type of response u which the circuit mil make, 
it is therefore only necessary to express u in the form 

and examine the bottom line of the fraction F{p). 


6*2. Resonance. 

A simple example is afforded by the phenomenon of Resonance, 
which occurs when an alternating stimulus is applied to a circuit, 
one of whose natural frequencies is equal to the frequency of the 
stimulus. An abnormally large response is apt to occur, as the 
reader may recall from steady-state a.c. theory. The following 
example illustrates what happens under transient conditions. 


6-3. Example. Resonance. 

A voltage ( F sin Qt ) . 1 is applied to an inductance L and condenser 
G in series, and = IjLC (Fig. 40). Find the current. 

The operational admittance of the circuit is 


1 Cp 

Lp + {llCpy “ LCp^+1 

^ ^ P 

Lp^ + Q^* 

Therefore the current is given by 


1 p 

Lp^+Q^ 


.{(F sin Dt). 1} 


I p VQp 


QV 


P 


L {p^+Q^y 
V 

JjJb 




Fig. 40. 


I 


i 


! 


I 


I 

i 

i' 


I 

i 
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The current therefore increases without limit, after the fashion 
shown in Fig. 41. 



6-4. Example. Resonance with Losses. 


In practice the current will not increase without limit, because 
losses ignored in § 6-3 limit its growth. Thus, if the inductance L 
has associated with it a resistance R (Fig. 42), the operational 
admittance becomes 


V 




Lp^ + { RjL) p + ( II LG) ’ Lp^ + 2ocp + ’ 

where a = a is supposed to be small compared with i?. 

The current is now 


I = 


QV 


P 




P 


.1 




L 2a \p^ + + 2ap + 


)■■■ 


F 

R 


I sin Qt--~ sin o)t \ , 

\ 0) J 


where The form of this is shown in Pig. 43. It tends 

to the final value 


F 

-^sini3^, 

Jti 


instead of to infinity. 
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The occurrence of resonance is signalized by the appearance of 
a squared factor of form + in the bottom line of the 
fraction F{p ) ; or, when losses are taken into account, by a product 
of form {p^ + {p^ + 2ap + Q^). 



6*5. Damping. 

The study of the p-expression F{p) in any problem enables us to 
determine a matter of great practical importance — the degree of 
damping which attends the oscillations of the circuit. Look back 
for a moment at the circuit of Fig. 9, whose different types of 
response are shown in Fig. 10. If the resistance B could be made 
zero, the voltage shock F would lead to a response v in the form of 
an oscillation persisting indefinitely. As B is made greater, the 
amplitude of the oscillation dies away more and more rapidly, 
until the ' critical state ’ is reached in which the oscillations dis- 
appear altogether (Fig. 10). If is made greater still, the re- 
sponse has no oscillatory tendency, but is represented by a com- 
bination of exponentials. This is typical of the behaviour of 
damped oscillatory circuits; a progressive change in the circuit 
resistances is associated with a progressive change in the degree 
of damping. 

It has already been observed that the clue to the behaviour of 
any particular circuit is to be found by examining the bottom line 
of the fraction F{p); this will have the form 

... + + 

and by finding the factors of this we can determine what kind of 
terms will occur in the solution, and whether oscillations are 
present or not. For every real factor (p + a), a term will 
occur; for every quadratic factor (p^+^p-i-y), where /?^<4y, 
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there will be terms like {B cos, (ot’^C mi (Dt). Thus we .can 
quickly determine the type of response which will be made by a 
circuit whose components are numerically known. 

But this is not enough. Often we desire to know, not only what 
the circuit will do, but how to make it do something else. It 
oscillates, and we want to stop it oscillating; or it refuses to oscil- 
late, and we w’-aiit to know how to encourage it to do so. Some of 
the values of the circuit elements are at our disposal; how must we 
dispose them to make the circuit do what we want? 

To answer this question, w’^e need a general criterion by which 
we can determine whether the factors of the bottom line of jF(^) 
are real or not. Writing this criterion in terms of the disposable 
circuit constants, we can find what relation must be satisfied by 
those constants in order that oscillations may be possible or im- 
possible. In Appendix V, where the properties of expressions like 

-h . . . 4 * Zp + m 

are briefly discussed, the reader will find these necessary criteria 
for expressions of the second, thii*d and fourth degree; the use of 
them will be illustrated by two examples, one of which is worked 
out very fully in order to illustrate certain features of this type of 
problem. 


6*6. Examples. Damping. 

Example 1 . A circuit consisting of L and G in parallel is shocked 
by voltage of unit function form, ap- 
plied through resistance R (Fig. 44). 

Find what condition R must satisfy, 
in order that the voltage v across L 
and C may oscillate, 

Theoperationalimpedanceof iand 
G in parallel is 

1 __ Lp Eig. 44. 

Gp + {llLpY LGp‘^'\~V 

Lp 



Therefore 


V 

V 


LGp^+l 

^LCp^+1 

Lp 

RLGp^ + Lp + R 



GET 


5 
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The condition for the factors of the bottom line not to be real is 
given in Appendix V ; in this instance, it leads to the result 


L^<4.R^LC, 


or 




which is the condition required. 


Example 2. A transformer is connected to a supply by way of 
a switch S ; but between them is an overhead line of inductance 1. 
Under certain conditions of imperfect operation of S, a large 
steep-fronted voltage shock V . 1 becomes applied to the terminals 
A, B; and by oscillation of the line inductance I with the trans- 
former capacitance G, the transformer experiences an oscillatory 
voltageof crest value 2 F. (The transformer inductance i', and the 
line capacitance c, exert a neghgible influence on this effect.) 

It is proposed to prevent this undesirable doubling by inserting 
a series reactor L damped by a resistance R (Fig. 45) . Clearly such 
an insertion cannot be effective if L is made too small. 



(1) How large must L be, if critical damping is to be possible? 

(2) If L can be made as large as this, what value should be 
given to i2? 

(3) If L cannot be made as large as this, what value of R will 
give maximum damping ? . 

( 1 ) By the usual method, 

V __ l/Op 

V~^{llGp) + lp-\-{RLpl[Lp^R])'^ 

Lp^R 

LlGp^ + RG{L 4- l)p^ + Lp + R^^^ ' ' 
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From Appendix V, the condition for over-critical damping is 

21 {Lie f - [EC{L + l)f + 4:{LlO) 

+ 4:[EC{L + l)f B ~ 1 S{LIG) EO{L + l)LE< 0, 

which reduces to 

+ {27 LWG ~-L%L + l)^G-l SL%L + 1) IG} E^ + {4:LH} < 0. 

Though at first sight this expression looks ugly, containing as it 
does so high a power as E^^ it is not so bad as it looks; for it only 
contains and so it may be described as being of the second 
degree in If L, as well as I and G, were specified, this inequality 
would enable us to determine whether any given value of R gave 
over-critical damping or not. 

Under what conditions can a quadratic expression become 
negative ? Think of such an expression as 

ax^-\-bx-\-c, 

where a is positive. The ‘roots’, or values of x which make this 
equal to zero, may be real (say/, g), or imaginary (say 
h—jh). If they are real, the expression has factors 

a{x-f){x-g), 

and this product is negative if, and only if, x lies between / and g. 
If they are imaginary, the factors are 

a[x — h—jh){x--h-\-jh), = — 

and this, being the sum of two squares, can never be negative. 
Thus a quadratic expression can only be negative if its roots are 
real and the variable lies between them. 

Our condition for over-critical damping can therefore never be 
satisfied unless the roots of the E'^ expression are real. The condi- 
tion for this is 

{21LWG -^L\L^lfG -I SL%L -h 1) Icf > 4{4(i + If C^} [4.1% 
or {27P ^ (L + Z)2 - 1 Sl{L -f l)Y > Ul{L + If, 

or - 24.m + 1%2LW - ^12LI^ > 0, 

or L{L-^lf>(^, 

or L>Sl. 

This shows how large L must be for critical damping to be 
possible. 


5-2 
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(2) If i is made equal to 8Z, the roots of the W expression are 
equal; there is then Just one value of R which gives critical 
damping. Substituting 8Z for L, we find that the expression 
reduces to 


12%C^R^-UmCR^-^4:Q%m = {S^CR^-2H)\ 


Critical damping occurs if this is zero; i.e. if 




?!1 


or R-= 1*54 

This answers question (2). 



(3) If i is made less than Si, critical damping cannot be 
attained. The expression for vjV can therefore have its bottom 
line factorized in the following way: 


V __ Lp + R ^ 

V^LlG{p + a){{p+/3f + o)^} 

For damping which is less strong than critical damping will 
permit damped oscillations (of the form e~^^GOS(j)t, e~^^sm.(jjt) to 
appear, and the corresponding factor in the ;p-expression is 
{{P + want to know how to make the damping as 
severe as possible — that is, we want /? to be as large as possible. 
Our first task, therefore, must be to find an expression for /? in 
terms of R, X, I and (7. 

Now 

LlG{p + a) {{p + 

was written for 

LlGp^ -{- R{L + 1) Gp^ + Lp + R. 

Therefore, equating the coefficients of p^, p and 1 on the two 
sides, we find 


and 
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From these three equations we can eliminate a and (j), so obtaining 
an equation for /?. The third of the equations gives 




jB 


whence, substituting in the second, we get rid of oi), and obtain 

cx.LlC~ W' 


In this may be substituted the value of a given by the first 
equation, namely (i?/Lo)~-2/?. The result, as we desired, is an 
equation involving /? only, which will be found to boil down to 
the following cubic: 

<W-2|(W+(|4)(¥)-i|-o. 

It must be constantly borne in mind that Lq, I and C are fixed, 
whereas we are able to select 
any value of R that we please. 

For each selected value of i?, ^ 
the real root of the cubic ^ 
equation determines a corre- ^ 
spending value of /? ; the curve 
relating the two will be some- 
thing like Fig. 46, for /? is zero 
when i? = 0 or CO, but takes 
positive values when R is in 
between. What we want to know is the value R^, which makes /? 
have its maximum value 

The cubic equation for (2/?) contains R to the second powder; 
in fact, it might be rewritten 

For any particular value of /?, say /?2 equation gives 

two values of R—R^ smd iJg. If the two values of R coincide, it 
can only mean that ^ = /?i; the condition for this is 




70 


RESO'KAIirCE, DAMPING AND STABILITY 

Therefore is the solution of the equation last written, which 
will be found to be 

4:U{L + l)C 

if we remember that = - y— - • . 

Substituting back in the i2-equation, we find 


L{L + 2l) 

2{L + l)^l{L + l)C' 


This is the value of R for maximum damping, if i < 8L 


It is interesting to notice that, if i= Sl,R^ = 


^ n 
27 J C 


or 1*48 


0 * 


Critical damping requires i? = 1*54 J Therefore, when L = U, 

critical damping is not quite the same thing as maximum damp- 
ing; though in fact the difference is trifling. 


6*7. Stability. 

All the problems so far considered have possessed a feature in 
common — namely, that we could make the maximum value of 
the response as small as we pleased, if we made the shock small 
enough. Usually the form of the response has been something 
like this: 


The terms on the right-hand side have a perfectly well-defined 
maximum value, say M; so u has the maximum value MU, 
which can be made as small as we please if U is sufficiently 
reduced. 

Only once have we seemed to encounter a different state of 
affairs, and that was in the first resonance example (§ 6*3). There, 
we found 


and it might appear that, no matter how small F might be, we 
should find that i attained a big value if we waited a very long 
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time. Blit this formula was obtained by neglecting the resistance 
which is inevitably present in every real circuit; and the suc- 
ceeding example showed that, when this is taken into account, 
the tendency for i to increase without limit is abolished. 

All this seems perfectly natural, until we recall that things are 
not always like this in other branches of physics. 

Think, for example, of a pendulum, consisting of a 
heavy ball on a light pivoted rod. It is theoretically 
possible for this pendulum to remain at rest in the 
upward position (Fig. 47); but let it be disturbed 
from, that position by the breadth of a hair, and it 
win swing right down. An infinitesimally small ' shock ’ 
produces here a large 'response’; for this reason, the 
upward vertical position is said to be ' unstable ’ . Can 
an analogous condition be found in an electrical 
circuit 1 

Let us -write down the equation of motion of the pendulum. 
We do not know the reaction in the rod, so we resolve the forces 
on the bob in a direction perpendicular to the rod, when this is 
inclined at an angle d\ the only force in this direction is the com- 
ponent of the gravitational force mg. Since 

Force in any direction = mass x acceleration in that direction, 
we have mg^ind = 



Until d gets too big, we may -write 6 instead of sin 6, and obtain 

If the initial value of 6 was 6^, we may -write the |)“equation , 


where _ =' gfL 

Hence, by the usual method, 


The term which indicates instability is for, no matter how 
small ^0 may be, it is possible to make quite big if we make t 

big enough. ; 
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In the electrical circuit problems so far treated, we have never 
encountered a factor like or cos cot; only cos cot, 

etc. Positive powers of e are the signal of instability. They never 
occur, however, in connection with simple networks of resistance, 
inductance and capacitance. Where they can occur is in circuits 
in which the value of one quantity is made to control the value of 
another — as, for instance, the grid voltage of a valve controls the 
anode current; the two examples given below will illustrate this. 
To detect the presence of instability, we need a criterion to de- 
termine whether any of the factors of a p-expression have the 
form (p — oc), or {p^—pp + 7), where a and ^ are positive; for such 
factors will lead to positive powers of c. For expressions up to the 
fourth degree, the reader will find these criteria in Appendix V. 

6*8. Examples. Stability. 

Example 1. Fig. 48 shows the circuit of a valve oscillator. If the 
{L, C) circuit is set in oscillation, an oscillatory voltage is applied 
to the grid of the valve by way of the mutual inductance M, This 
causes pulsations in the anode current of the valve, which can — 
under certain conditions — overcome the natural tendency of the 
circuit oscillations to die away. It is required to find these 
conditions. 


(volts) \ i 
^ A B 

Fig. 48. Pig. 49. 

The characteristics of a typical valve are shown in Fig. 49. 
The anode current is dependent jointly upon the anode voltage 
and the grid voltage and over a considerable range AB the 
‘ curves ’ are straight lines, so that we may write 

where a and 6 are positive. 
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Suppose tliat the circuit is miming steadily without oscillation, 
drawing a direct current from the n.c. supply . The quanti- 
ties and must then be related by the equations 

and •S^ao + eao = 'fo- 

Now suppose that something happens to change the current 
in (Jy, B) to the new Talue ^^0 + ^ voltage across C to 

Bi^o + v, We can obtain a differential equation for v and i by 
virtue of the fact that 


Voltage across 0 = Voltage across L and B; 


that is 

Riao + V = + R{'i'aO + ^)> 

or 

y. d/i , 

V == L-7- + 

at 

We can obtain a second relation from the characteristics of 
the valve. The steady anode current, anode voltage and grid 
voltage now change to variable quantities e^, related to the 
original quantities as follows : 


= Current in i-f Current in C 


. ^ dv 


~ Voltage across C 


--d% 

But 

\i — ^*a0 ~ "h bCgQ + C. 

Therefore 

C-y+% = -av + bM-j-. 
at dt 


We now have two differential equations linking v and i with t 
The corresponding ^-equations are 

v — (Lp + R) i = terms depending on initial disturbance, 

{Cp + a)v~{bMp — l)i — „ „ „ 
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If we are only interested in finding the conditions for continued 
oscillation, we need not worry about exactly what terms occupy 
the right-hand side; for these will only appear in the top line of 
the _p-expression for v or i, whereas the presence of oscillations is 
made evident by the bottom line. So, eliminating v (for instance), 
we get 

{{Gp + a) {Lp -f i^) — bMp + 1} i = initial terms. 

The conditions for sustained oscillation can be found by ex- 
amining the p-expression in curly brackets, which is 

LGp^ + {RC^aL-hM)p + {l-\-aR). 

We need only refer to the criterion of stabilit}?^ given in Appen- 
dix V ; in this case the 'unstable ’ condition is the condition to be 
desired. Having observed that the preliminary condition, 
LG > 0, is satisfied, we write down the conditions for the system 
to be stable, namely 

(1) RC -\-aL — bM >0, (2) l + ai?>0, 

and observe that there will be instability if one of these conditions 
•is broken. But condition (2) cannot be broken; so instability is 
only possible if 

RG-^ctL — < 0, 

This is not quite enough; for an oscillator is of little use unless it 
oscillates. The condition for there to be any oscillations at all must 
also be satisfied, and this is 

{RC^aL- bMf < 4LG( 1 -f aR), 

Failing satisfaction of this second condition, the quantities i 
would simply increase like e^\ instead of cos o)t, etc. 

If both these conditions are satisfied, the current will be found 
to be of the form 

i — A^-e‘^\BQOB(i)t+Gmi(j)t), 

which might seem to betoken an oscillation of ever-increasing 
amplitude. What actually happens, of course, is that the ampli- 
tude increases only until the curved parts of the valve charac- 
teristics come into the picture; the assumptions of this discussion 
then cease to be true, and instead of continuing to increase the 
oscillation settles down to a steady maximum value. 
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Example 2. The voltage of a d.c. generator is intended to be 
held against changes of speed, by a regulator of the form shown 
in Fig. 50. A coil, energized by the generated voltage, attracts an 
armature against a constant counter-pull Tq\ the armature is 
coupled to the moving contact of 
a variable resistance in the field 
circuit. A momentary increase 

in generated voltage, due to an --- 

increase of speed, should there- 
fore be automatically corrected 
by an increase in the resistance of 
the field circuit. 'Armature reac- 
tion ' in the generator is neglected. 

The movements of the regulator Fig. so. 

are damped by an oil dashpot (not shown). It is required to find 
how much damping must be provided, in order to prevent the 
system from oscillating persistently. 

It will be observed that this system can only be in equilibrium 
when the generated voltage has a particular value Vq, namely that 
voltage for which the solenoid pulls the armature with a force . 
If the speed has its normal value we will suppose that a field 
current suffices to generate this voltage. Then the resistance 
in the field circuit must be given by 

Vq = 

Now suppose that the speed suddenly changes to -f 9^, where 
njuQ is small, and remains constant at this new value. Everything 
else starts to change in sympathy. At a subsequent instant, let the 
generated voltage be + and the field current i^ + i; further- 
more, let the regulator armature move a distance x from its 
initial equilibrium position, increasing the field resistance by a 
proportional amount which we will call ax. 

Our first equation will come from the mechanical behaviour of 
the regulator armature. If the extra voltage v causes an un- 
balanced force bv, this equation is 

d^x jdx , 

where m is the mass of the armature, and J; is the damping 
coefficient associated with the dashpot. 



Fig. 50. 
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Secondly, we have the relation between voltage and current 
in the field circuit; this is 

d'b 

where L is the self-inductance of the field. Neglecting the product 
axi (for both x and i are small), and remembering that = Vq, 
we may rewrite this equation in the simpler form 


_ U/D ^ . 

4 * + aiQX = V. 



The third equation involves the saturation curve of the 
generator. This will have the form shown in Fig. 51. At constant 
speed tIq, the increase v in voltage is related to the increase i in 
field current by the equation 

V = ci, 

where c is the gradient of the tangent CG' to the saturation 
curve; but the speed variation changes this to 

(vo+?;) = (t)o + a)(l+|^j, 



neglecting the small quantity (wci/Wfl)- 

These three equations (two differential, one ordinary), de- 
termine the behaviour of the three unknowns v, i, x, in terms of 
the known speed change ». 
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The corresponding ;p-eqiiations are 

{mp'^-{-kp)x — bv = 0, 
ai^x — -f [Lp i = 0, 

v — ci = {nv^TiQ). 

As in Example 1, there is no need to know what terms occupy 
the right-hand side; bnt since v is the only quantity which can 
change suddenly with a sudden speed change n, the right-hand 
sides must actually be as shown. 

Eliminating x and we find . 

7W 

{(Lp -t (mp^ + hp) -h ahcijj} v = {Lp + 5^) + ^P) — ~ ? 

Uq 

w^here Ri = Bq~ c. Referring to Fig. 51 again, we observe that Rq 
is the gradient of the line RR'; so R^ is the difference between the 
gradients of RR' and CC'; it is positive. 

The stability of the system therefore depends upon the 
properties of 

{Lp -f J?i) {mp^ + kp) + abciQ, 

or {Lm)p^ + {Rpn-{- Lk)p^ + {R-Jc)p + {abciQ), 

Referring to Appendix V, we find that the conditions for this 
system to be stable are as follows (the p-expression has been 
written, as it must be, with the coefficient of p^ positive) : 

(1) Bpn + Lk>0. (2) abci^>^. 

These two are automatically satisfied, for a, b and c are all 
positive. 

(3) R-Jc{Rpra -}- Lk) > LmabciQ. 

This last equation gives the 'safe’ values of k. It may be re- 
written 

R^Lk^ 4 - Rlmk ~ LmabciQ > 0 , 

and the only positive values of k which enable this to be satisfied 
are those for which 

, -R^m-i-^lRfm^ + 4R^L^7mbciQ 

■ 

This is the required criterion of stability. 
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. But absence of sustained oscillations does not in itself mean 
that a control system is satisfactory. Usually the engineer re- 
quires that the controlled quantity shall attain its final value in 
a critically damped manner; or perhaps he permits a little oscil- 
lation, because this allows the, quantity to get near to its final 
value more quickly (Kg. 52). A convenient supposition, which 
gives about the right degree of damping, is that each oscillatory 
term cos (i)t or siaojt is damped by a factor condition for 

this is also given in Appendix V. Applied to the present problem, 
this condition leads to the criterion: 

{Lmf d? - (R^m + Lkf {RJkf + 2(Lm) {R^hf 

-+■ 2d (R^m -f Lhf — 4(im) {R^m + Lk) {R-Jc) c? = 0, 

where d = abciQ. Rewriting the criterion as an equation for k, 
we obtain 

- Rim^ + 2mB - {R{7n^ - 2R^L^md) B 

-f 2R\Lm^dk 4- + 2R\m^d) = 0. 



This equation of the fourth degree has a real positive root, which 
will be found to be somewhat larger than the value of k w^'Mch 
ensures stability; large enough, in fact, to ensure the degree of 
damping which is required. 

If the condition for critical damping is written down, it leads 
to the rather similar equation 

-- Rim^ + {4.m 4- 2R\Lm) F - + 6R^L^md) B 

-^RllMiHk+{21L^mW^mimH^ 

and the positive real root of this equation will be larger still, for 
this is the stiffest of the three criteria which we have discussed. 


CHAPTER VII 


SOME PRACTICAL EXAMPLES 

7*1. The reader of the preceding six chapters is now the 
possessor of a set of tools, new and shining, for prising open the 
mysteries of electrical circuits. But the circuits which we have 
discussed have been chosen so that each displayed chiefly the 
prowess of the newest tool; we have not yet had to lay the circuit 
on the bench while we scratched our heads dubiously before the 
tool-box, wondering which tool to select. In the examples of this 
chapter we shall pass through this final discipline. They are 
chosen from several branches of electrical engineering, and in 
each we shall have to consider what is the best of the methods 
available to us. 

7-2. Ideal and Non-Ideal Interruption of Alternating 
Current. 

Two taps on a transformer winding are connected by contactors 
to the extremities of a reactor, which is used in on-load tap- 
changing (Eig. 53). The reactor has inductance 100 millihemies, 
and equivalent shunt capacitance 1000 micromicrofarads; its 


Transformer 




losses are equivalent to a parallel resistance of 50,000 ohms. The 
transformer reactance is negligible by comparison. Find the 
wave-form of voltage which appears across one of the contactors 
when it interrupts a current of 30 amps., 50 cycles, 

(a) exactly at current zero; 

(5) before current zero, when the instantaneous current value 
is 0*2 amp. 

The relevant portions of the circuit are set out more clearly in 
Fig., '54. 
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The method of Impedance operators is iisnally the quickest, 
unless the allowance for initial conditions is difficult. Here, the 
voltage which we want is zero until the instant of switch opera- 
tion; so the value afterwards is the same as would be produced by 
injecting across the switch terminals into a dead circuit a current 
equal and opposite to that through the closed switch (Fig. 55). 

Intelligent simplification should never be despised. We observe 
that the frequency of the (L, G) 

circuit is 15,900 cycles, whereas [ ~T | ~1 

the supply frequency is 50 /stnntl i c ^3 

cycles. For many cycles of the , "T > S 

(L, G) oscillation, therefore, ! f r 

the injected current will not 
diverge perceptibly from the 

straight line IQt (when interruption is at current zero), or IQt — Iq 
(for pre-zero interruption). This simphfies the working con- 
siderably, without any loss of accuracy so far as the interesting 
part of the transient period is concerned. 

The operational impedance of the circuit is 


Cp + illB) + {llLpy 


RLp 

RLGp^ 4- Lp + i? ’ 


If we work out problem (6), we can treat {a) as a special case. The 
voltage produced when (IQt — lQ ) . 1 is injected into this circuit is 
given by 

RLGp^^-Lp-^R ^ 

_ RLp ilQ \ 

” RLo^^Lp + i? ‘ -^7 ‘ ^ 

_ -- JqP + IQ 

Gp^ + {llRG)p + {llLG) 

Since the next step depends on whether the factors of the 
bottom line are real or complex, it is best now to substitute for 
i?, L and G, Measuring time in microseconds, we obtain 

V = 1000 - + 1 

^3^ + 0-02p + 0-01 

= 1000^ — + 

CP + 0-01)2+ (0-0995)2 ’ 
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but in most work of this type we are not interested?in the dif- 
ference between 0-0995 and 0- 1 , so we write 

V = 1000; — + 

(p -t- 0-01)2 -I- (0-l)2-'‘- 

Turning now to the table of standard forms, we find that the 
interpretation is 

V = 1000 1 — ^ g-o-ou gjjj 

+ + )• 

In problem (a), = 0, 1 = 30^2, i3 = 3-14 x 10~^ (when time 

is measured in microseconds). Therefore 

V = 940V2{l-e-«-®“[cos(0-l<)-F^sin(0-l«)]}. 



Fig- 56. Fig. 57. 


The form of this voltage is shown in Tig. 56. The oscillation is 
dying down to the steady value 940^^2, which is the crest 50 cycle 
voltage across the open contacts. This is natural, since from the 
power frequency viewpoint the circuit is practically pure in- 
ductance. 

In problem (6) there is superimposed upon this the contribution 
of the If) term, which is 

— 2000e“"'®^^sm (0-l<). 

The resultant shape is shown in Mg. 57. It makes a lot of dif- 
ference whether the current is extinguished at 0-2 amp. or O-Oamp. 

CBT 6 : 
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■ 7-3. Recovery 'Voltage upon Extinction of an Earth Fault in 
a Supply System Earthed through a Petersen CoiL 

A three-phase star-connected transformer supplies a trans- 
mission system, each of whose lines has capacitance K to the 
other lines, and (7/3 to earth. The transformer neutral is earthed 
through an inductance L (known as a Petersen coil), whose value 
is such that l/27r^/LC is nearly equal to the supply frequency. 
The losses in the coil and system can be represented by a resistance 
R in parallel with the coil (Fig. 58). Compared with L, the trans- 
former reactance is negligible. 


, Transformer 
(Primary not shown) 


\ kJ 

-X- : 

rX 



rT~v ^ 

i ^ / . 

± 

XT 



ii. F'k 
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9.^ 

TTTTTTf 

Yrrf 

. 1 


Fig. 58. 


At the point F an earth fault has become estabhshed, for a 
sufficient time for the current in the fault to have attained a 
steady (alternating) value. This current becomes extinguished at 
a zero point; we desire to know the wave-form of voltage appear- 
ing across the fault. 

This problem is made to look much simpler if we replace it by 
the equivalent one in which a current equal and opposite to the 
fault current is injected into the circuit at the fault vnth all 
constant voltage generators regarded as having zero impedance. For 
the transformer here plays the part of a constant voltage gene- 
rator; if it is assigned zero impedance, the capacitances K are 
short-circuited and disappear from the problem, while the 
capacitances CjZ are connected in parallel. Thus the problem 
can be represented in exactly the form shown in Fig. 55. 

This makes the problem look very like the last one, with the 
significant difference that the natural frequency of the circuit is 
now nearly the same as the supply frequency. We could start by 
working out the value of the steady fault current; but we save 
ourselves this labour by calling it I sini7i5, and resolving to give / 
such a value that the ultimate voltage equals the transformer 
phase voltage— say F cos(i3^-j-0). 
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The operational impedance of the circuit is 


RLp 

BLCp^ + Lp + B' 

Therefore the voltage produced by injecting / sini3i . 1 is given by 

^ ~BLCp^ + ip ^ • {-f sin . 1} 

IQ ^ 

where 2a = l/iJC?, a^ + = 1 jLC. 

Splitting up the p-expression into partial fractions, we suppose 
_ Ap^ + Bp Ep^+Fp 

[(p + a)2 + «^]ljp2 + i32]~(p + a)2 + (y2+ + 

then, multiplying both sides by [{p + a.f + [^2 + ^ 2 ]^ obtain 

A = —{(x^ + (a^—Q^)ID, where X) = (a^+w®— + 

B = —2ci{a^ + (t)^)/D; 

E = {a^ + o}^-D^)ID; 

F = 2aQ^ID-, 

G = 0. 

Substituting these values, and performing the p-operations with 
the aid (if necessary) of Appendix I, we find 

V = |(a^ + (x)^—D^) (cos Qt - 6““' cos mt) 

+ 2ai3 sin - - (a^ + + Q^) e~“' sin at ■. 

We have still to find the value of / in terms of V. When ^ 
becomes very large, v tends towards the value 


IQ 


{{a^ + (o^-Q^) cos Qt + 2aQ sin Qt), 


IQ 

G^D 


cos {Qt+ 6), 


where cos 6 = {oc^ + a^— Q^)I^JD, sin 0 = — 2aQl^jD ; 
the bottom line has of course been adjusted to make 
cos^0 + sm®0 =1. 
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Therefore IQjC^D must equal F ; and we may rewrite v in the 
form 


'y = Y 


I (a^ + Q^) (cos — e”«*cos a>^) + sin Qt (a^ + H- Q^) m 


V(a^ + 

A more seemly form is 

'D = F{cos {Qt + 0) ~ Ke~^^ cos {o)t + 
where, it will he found, 

cos^ == (jt){a^-}-(i}'^—Q^)l'jD{a^ + o)^), 
sin^ = —oc{(X^ + o)^ + Q^)l^jD{a^-\-o)^). 




To gain an idea of the form of this voltage, we must bear in mind 
that a/w is quite small — ^not. greater, probably, than 0*04. For a 
preliminary survey, therefore, we put a = 0, and obtain 

= F(cosi3^-coscL)^) 

. (Q + o) \ . (Q-a)\ 

Since w is nearly equal to Q, j represents quite a low 

frequency compared with - j; and the appearance of the 

wave is therefore that of a wave of frequency ^ x j , whose 

amplitude varies at the frequency ^ x V The mavirrinTn 

27r \ 2 / 

value is nearly 2F (Mg. 59). The actual wave wil start as in 
Fig. 59, but the pulsations will gradually fade out into a steady 
50 cycle voltage (Fig. 60). 


150 

YDS. 


300 

YDS. 
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7-4. Modification of a Lightning Surge Voltage by the Cable 

system of a Power Station. 

trari 1 ^ cables, respectively 150 and 300 

yards long, connected at A and C to • 

overhead lines ; BD is a cable 10 miles ® 3:~ ^ B 

in length. 

A hghtning stroke to the overhead 
line near the point A causes a high 
voltage of approximately ‘unit fimc- 
tion shape to be apphed to the system 

at A. It is required to find the voltage 

on the busbars B, under the following assumptions: 

r, behave approximately as though 

half their capacitance were concentrated at each end, with their 
inductance between; actuaUy, of course, these are distributed, 
a little bit to each bit of cable. 


1 


Pig. 61. 


(2) That the cable AB has a total series inductance of 30 
microhenries, and a total capacitance to earth of ^ microfarad- 
and that the other cables have the same inductance L and 
capacitance 0 per unit length as AB. 

(3) That the cable BD, which acts more or less as a continuous 
absorber of energy, may be treated as equivalent to a resistance 
of ohmic value \j{LIC): (With the values of L and C already 
given, this comes to 30 ohms.) 

(4) That the overhead fine at C exerts a negligible effect on the 
behaviour of the cable, i.e. C behaves like an open circuit. 



3QaH 



— — 

•Q5/^F ^ 


Kg. - 

60>ttH -033^? 

63. 



Under these assumptions, we shall be replacing the cable 
system by the network shown in Kg. 62; this is redrawn in a 
clearer form in Kg. 63, where the capacitance between the 
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point A and .earth, has been omitted, since it is in parallel 
with the whole system, and makes no difference to the voltage 
distribution. 

The operational impedance Z between the point B and 
earth is given by 

if we work in microhenries and microfarads; which reduces to 

3 ^^ + 2 ^^ + 2 * 5 ^ + 1 

30(2jp2+i) • 

Therefore the voltage v at the point B, worked out just as though 
this were an ordinary a.c. circuit, is given by 

T7 ^ 1 

^BOp + Z'^ 

= F .1. 

To interpret this expression, we need to find the factors of the 
bottom line. The method is given in Appendix V. 

We write 

3 ^^ + 2 ^® 4 - 4 - 5 ^^ + 1 

in the form 

3{(i5 ® + ip +/)2 - {gp + h)% 

and proceed to discover suitable values of/, g and A. Equating 
corresponding terms on the two sides, we find 

2/-^^ = ft. 

The first and last of these equations give g^ and in terms off; 
the second shows that ZQg%^ = {2f—l)^. We can therefore 
eliminate g and h and obtain an equation for /, which turns out 
: tobe 
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Tlie first equation also shows that g will only be real if/ is greater 
than If; of the three roots of the cubic equation for/, only one 
satisfies this condition; the equation is solved as an example in 
Appendix V, and the root is there found to be 0*709. We can then 
find at once that g = 0*169, h = 0*411, so 

3|)^ + 2p^ + 4*5p^+p-l- 1 

= 3{(p2 + + 0*709)2 (0*169^ + 0*41 1)2} 

= 3(p2 4- 0*502p + l*120)(p2 4.0*164^ + 0*298). 

These two quadratic factors cannot be split into real factors of the 
first degree. 

We now seek partial fractions, in the foEowing form: 

2p2+ 1 

3(p2 + 0*502p +“1*120) (p2 + 0*164p + 0*298) 

_ Ap^^-Bp (7p2 + i)^ 

p2 o*502p + 1* 120 p2 4. Q. 104^ ^ 0*298 * 

Multiplying both sides by 

3(p2 + 0*502p + 1* 120) (p2 + 0*164p + 0*298), 

and equating corresponding terms, we find 

A = --0*533, 5 = -0*064, 0 = -0*467, D = -0*280, = 1. 

(The last named may be more simply obtained by putting^ = 0.) 
Therefore 


0-5d3p^ + 0-064p 

0-467jp2 + o-280p ) 

(j) + 0-251)2 +(1-028)2 

(p + 0-082)2 + (0-540)2/ 


= ys^l - e-0‘25i*5(0*533 cos 1*028^-0*068 sin l*028s5) 

- e-<>‘^s2^(o.407 eos 0*540i -f 0*448 sin 0*540^)}, 

which is perhaps more conveniently written 

= F{l-0*537e-o-25Ucos (1*028^ + 0*127) 

- 0*647 cos (0*540^ - 0*764)}. 

The form of this voltage is shown in Fig. 64 (full curve). Its 
maximum value is 1*35 times the applied voltage. 

Now suppose that the applied voltage, instead of being of 
'unit function ’ shape, rises steeply to its maximum and then 
falls to half that value in 20 microseconds; what happens to the 
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voltage at B ? This solution may be deduced from the other with 
the help of DuhameFs integral; The applied voltage, instead of 
being F. 1, will now be x 1. 

If we apply the Duhamel method in the obvious "way, we shall 
reach a solution something like the^ solution already found. 
Plotting such solutions is tedious, as the reader will discover if he 
tries it. Since we have already plotted this solution in Fig. 64, 
why not deduce the other by numerical integration — an easy 
arithmetical job? 




The essential of the Duhamel method is to represent the shock 
on the circuit by a series of "unit functions’. In this instance, it 
will suffice if we apply the successive "unit functions ’ at intervals 
of 2 microseconds; but, in order to make the built-up curve zigzag 
on both sides of the true one (see Fig. 65), we apply the kep, which 
lowers the curve from its O-microsecond to its 2-microsecond 
value, at time 1 microsecond. Subsequent steps then occur at 
3, 5, 7, ... microseconds. Omitting the constant factor F, the 
staircase which we use to represent x 1 is therefore as 

follows: 


At time 0 microseconds, 


+ 1 . 


>5 

J5 


1 

3 

5 

7 


55 


>5 

55 


_(l_e-o-o7) =-0-0676. 


_(e-0'07_e-o-i4) = 

-0-063o. 

_(g-0-14_g-0-21) ^ 

-O-OSSg. 

_(e-o- 2 i_e-o- 28 ) = 

-0-054g. 

_(e-0-28_g-0-S5J _ 

-O-OSlj,. 
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Tlie response of the system to this staircase is as follows : 





(Time ttsec.) 



Component 







I 

2 

3 

4 

5 

6 


1 

0*252 

0*666 

0*921 

1*007 

1*071 

1*191 

- 0*0676 

0 

- 0*017 

- 0*045 

- 0*062 

- 0*068 

- 0*072 

- 0 * 063 o 

— 

— 

0 

- 0*016 

- 0*042 

- 0*058 

- 0*0583 

— 

— 

— 

— 

0 

- 0*015 

- 0*0548 
- 0 * 051 i 


— 

— 



— 

Total 

0-252 

0*649 

0*876 

0*929 

0-961 

1*046 

Component 

7 

8 

9 

10 

11 


1 

1*315 

1*342 

1*236 

1*048 

1*100 


- 0*0676 

- 0*081 

- 0*089 

- 0*091 

- 0*084 

- 0*071 


- 0 * 063 o 

- 0*063 

- 0*067 

- 0*075 

- 0*083 

- 0*085 


- 0*0588 

- 0*039 

- 0*054 

- 0*059 

- 0*063 

- 0*070 


- 0*0548 

0 

- 0*014 

- 0*036 

- 0*050 

- 0*055 


- 0 * 051 i 

— 

— 

0 

- 0*013 

- 0*034 


Total 

1*132 

1*118 

0*975 

0*755 

0*785 



The first component is read off from the full curve of Fig. 64. 
The others are simply obtained by multiplying the figures in the 
top line by the appropriate factor, and also shifting them to the 
right by the required time-interval. The whole process only re- 
quires a few minutes of slide-rule work, once the response to unit 
function has been calculated; the total, giving the response to 
F X 1, is plotted as a dotted curve in Fig. 64, . 

We have been using DuhameFs integral in the form 

= C7(0)%(^)+J U'{r)u^{t-7)dr, 

or U{Q)uS)^i^\r)ST.uS-r). 

0 

Z7(0)^i( 35) is the first component; U'(t)St is the height of the 
^step’ — which, however, we have located at the instant (r — 1), 
so that instead of Ui{t — r) we have used — r -{- 1 ) . 

7-5. Behaviour of a Six-Phase Rectifier Circuit. 

A six-phase rectifier equipment of negligible internal impedance 
supplies the terminals A, B of a smoothing circuit {L, C), which 
is loaded with a resistive load B (Fig. 66). The voltage applied to 
A, B is the envelope of a series of 50-cyele sine-waves, of crest 
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value F, displaced by 60® (Fig. 67); but the rectifier can only pass 
current in tbe direction shown by the arrow, so that if the current 
tends to reverse it becomes extinguished. The terminals are then 
effectively disconnected from the voltage supply until the voltage 
across them falls lower than the instantaneous value of the 
supply voltage, when the arc re-strikes. 


For very light loads, the voltage across G and M will remain 
almost constant at the value F. For heavy loads, the arc will 
remain in being throughout the cycle, and the mean voltage v,^ 
across C and R will equal the mean value of the rectified voltage, 
which is 3F/7r. We shall therefore ^ . p 

expect to be related to the mean 

by some such curve 1 


load current i 
as is shown in Fig. 68. It is required ^ j 
to find, for the values of L and C I 

shown in Fig. 66, the point P at | 

which the voltage starts to rise; ; 

that is , t o find the load at which the C 

rectifier current just touches zero. 

It is perfectly possible to obtain an exact formula for the wave- 
form of rectifier current, by using the method set forth in § 4*7, 
Example 2. But unless the reader has a taste for several hours of 
algebra, he had better start with a bit of intelligent simpHfication, 
as follows. 

The voltage wave shown in Fig. 67 oscillates between the values 
F and 0*866F , at a frequency of 300 cycles. For a rough approxi- 
ruation, replace the wave by a sinusoidal ripple of the same 
maximum and minimum, and the same frequency (Fig. 69). The 
formula for this wave will be 

F(0*933-i- 0*067 sin 6i30, 
where = 27r x 50 

= 314, if time is measured in seconds. 


where 
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At 300 cycles, L has impedance 3014 ohms; G has 26 ohms. 
Therefore, whatever the value of R, the 300-cycle impedance of 
the circuit must lie between 3014 and 2988 ohms. Calling it 3000, 
we see that the peak alternating current through L is 


(0-067F/3000) amps. 

The direct current depends only 
on R, and is 



(0-933F/i2) amps. 

If the direct current is just large Pig gg 

enough to prevent the current from 

tending to reverse, these must be equal; from which we find 


R — 42,000 ohms (approximately). 


The critical damping resistance of the (i, G) circuit is 


=141 ohms. 

Rough though this method is, it makes it quite clear that the 
interesting value of R is far larger than the critical resistance; 
so much larger, that it exerts only the shghtest influence on the 
oscillations of the {L, G) circuit, and is reaUy effective only as a 
drain for direct current. This suggests the following method of 
attack : 

(1) Find the current through i when the wave-form of Mg. 67 
is applied to the circuit of Mg. 66, with R open-circuited. 

(2) Mnd the greatest negative value of this current, and hence 
ascertain what value must be assigned to R so that the direct 
current just cancels this negative a.c. peak current. (For this 
purpose, the voltage across R may be assumed to have the mean 
value of the voltage form in Mg. 67; namely, SF/tt.) 

In calculating the current when R is open-circuited, it would 
not be wise to use the method of § 4-7, Example 2; for, in the 
absence of damping, the effect of each voltage cycle lasts for ever, 
so that the cycle which happened an hour ago may have a residual 
effect as large as the cycle which has just been' completed. We shall 
adopt a method ^^hich is less elegant, but more safe. 
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Take as variables the current i in L, and the voltage v on C 
(Fig. 70), The applied voltage may be taken as F s,m{Dt-h 7 r/Z), 
during the interval from ^ = 0 to ^ = n/W, 

At the beginning of the interval, let i = L, 

and^; = ^;„. 1 


Writing down the 7 }-equations by the , n, 

i-i, ^ e /LX ^ lA • ¥svn{at+f) 

method of § 5*5 (6), we obtain 


Lp{i-iQ) + v = VBm(Qt-\~7T/Z), 
Gp{v — Vq) = i. 


Fig. 70. 


Eliminating v between these, 

{LCp^+l)i = LGpHQ—GpvQ+CpVBiii{Qt-^7rl3), 

or (p2 ^ 4- ^ 7 sin {Qt + 7r/3), 

Jj Jb 

where wl^—ljLC, 

Thus, writing V sin(i3^-h7r/3) in its oiierational form, we obtain 
the operational form for i\ namely, 

• f % P V ^Zf^ + Qp^ ] 

Writing the last term in partial fractions as usual, and inter- 
preting, we find 

^ cos (i)t — 

Lo) 

V 

+ ” COS Ot) - o>^3 sin o)t + Q^Z sin Qt}, 

'< We now have to find and Vq by writing down the conditions 
that i and v may have the same values when ^ = 0 and t = 
(where = tt/ZQ); for the circuit must be assumed to have 
attained a steady state. These conditions are 

i = Iq when t = 


idt = 0: 


SOME PRACTICAL EXAMPLES 


93 

for the integral of i over any interval of time is equal to the charge 
acquired in that time by the condenser. In this way, we find 

1 — cos sin (ji)U 

JUO) 

V 

== 2L { Q^ cos -0)^3 sin 

V 

and t’o sin ^ ( 1 — cos dJtj) 

V 

^ 2L(Q^-d)^) ^ ^^i}- 

Solving these equations, we obtain 

i’Q = 0 , 

o)V Q^m(j!)t^ — o)^3{l — cos(i)t^) 

2(Q^ — (ji>^) (l-~cos(t>q) 

_ (j)V Q cos {(i)t J 2) -- o)^j3 sin 2) 

2(Q^ — o)^) sm{o)til2) 

Hence, after a little trigonometr^’^, we find 

i3F f . (o)ti \ (ot^ ( 7r\) 

A simpler form may be obtained if we measure time from the 
instant when the applied wave reaches its crest. Using electrical 
angles instead of time-units, let 

Qt — 71/6 = 6^ 

and (jt)IQ = A. 

The expression now becomes , 

V i, sinA^ ] 

^ ~ 2sin(A7r/6))’ 

where 6^ varies from —7r/6 to + 77/6. 


94 


SOME PRACTICAL EXAMPLES 


With the particular component values given in Fig. 66, 
w = 177, so A = 0-564. Therefore 

i l-7183sin(0-664(9)}, 


and the form of this current is plotted in Fig. 71. The negative 
crest occurs when 6 is about — 17*^, audits value is — 0*00927 F / LQ. 

The mean direct current through L is ZV/ttM, If this is just 
enough to nullify the negative a.c. crest, we must have 

R = 51,800 ohms. 

It would of course be possible to allow for the fact that the voltage 

across R is really v, the voltage 

occurring in our equations; but .oi— 

it is doubtful whether such re- 

finements can be justified in this / \ 

approximate method. Anyway, 3% \ q 

the instant when i attains its 
crest must be quite close to the 
zero of the oscillation on so the 
instantaneous voltage across R 
will be very near to the mean 
value of V, which must be SV/n. 



Fig. 71. 


7*6. Surge Propagation in Alternator Windings. 

One phase of an alternator winding consists of 24 coils connected 
in series, each with self-inductance L and each having a capaci- 
tance C to the earthed core. The capacitance is really distributed 


Y 


Fig. 72. 


IsrCOIL- 


-2nd COIL- 


hva Va I I 


-24th COIL- 




EARTH 




T 


along the conductors of the coil, but with sufficient accuracy it 
may be considered to be concentrated at the centre of each coil, 
with half the inductance on either side of it (Fig. 72). If a surge 
voltage of 'unit function’ form is appMed between terminal A 
and earth, what is the form of the voltage to earth at B, the 
junction of the first and second coils? 
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It is to be hoped that the reader will not attempt to write down 
the complete ^-equation for this circuit, since it is of the 48th 
degree. It is useless to attempt a problem like this unless some 
simpler method can be found; for instance, the following. 

The winding between the points A and ‘Earth’ will have an 
operational impedance, Z. The last 23 coils, to the right of the 
point B, will have practically the same impedance; for what is 
one coil among so many ? The condition for the two impedances 
to be equal is 


L'p 


+ - 






Lp 


+ Z . 


and this is an equation for Z, whose solution is 
Z-> 


/i LC A 

(i+Ty-)' 


(The reader will see that what w'-e have really done is to assume 
that the number of coils is infinite.) 


The circuit has now been simplified 


V 2 


u 


to that shown in Fig. 73. The current 

V V \ 

in AM is 7 =.!; in MB it is and 

Zj Z \ 

the difference of these is the current 

taken by the condenser C, Fig. 73. 


5 


L/2 


T 

V 

i 


= GpVm. 

v^{l~ZCp + iLCp^)VA, 

This bears no resemblance to any ‘standard form’ which we 
have discussed so far. Nevertheless we can get quite a long way 
by simply applying our known rules of interpretation — that is, 
by expanding in powers of 1/j^ and obtaining the solution in the 
form of an infinite series. 


Therefore 

or 
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Substituting for Z, we find 

« = F(1 + - 4lG^{1 + W^)) ■ i ■ 

Writing l/w* for LG, as we have usually done, 


« = F 1 + 






By taking out the factor p^j2(d^, we may transform the bracket 
to a form which contains Ijp instead of jj; thus, 


Fp2 

2w2 


1 + 


p2 


1 +- 


4=0)^ 

p2 


. 1 . 


Now use the Binomial Theorem to expand the square root in 
powers of (4w^/p^) — as explained in Appendix II. We then obtain 


Fp2( 2 c<j 2 r 14m^ 1.1 Uo)^ 1.1.3/4w2\» 1\ 

^ ~ ■^23.3! •••Jj' 


1 , 


Fpa 

26j3 


1 1.3ptoY i_z.5 l2o}^Y 

2 ! \ j 3 ! \ j 4 ! ^ 


p- 


. 1 . 


The general term of this series, including the factor p^l2o}^, is 


•1)’ 


1.3.5...(2ra-3) 
n ! 


\p2 j 


If we multiply the top and bottom Mnes by 2.4.6 (2n — 2), 

which is 2™“^ x (» — 1) !, we obtain the simpler form 

(2n-2)! 

^ ’ ^n-l)!w!\pV ■ 


Therefore u = F 


(»- 

2! "(Y 4! <u* 


[+■; 


6! w® 


1 ! 2 ! 2 ! 3 ! p* 3 ! 4 ! p6 

|(wi)3 {(atY ^ (cotY {(otY 
\1T2! ' 


.1 




2!3! 3!4! 4!6! 


+ .. 


Here we have obtained v in the form of a series of quite a simple 
type; the amount which we can discover about the behaviour of 
V is now only limited by our patience in working out the terms of 
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the series. Fig. 74 shows the behaviour of v up to the time given 
by = 5. It shows that the voltage at B will reach a crest value 
which is about 13 per cent in excess of the applied voltage. 

This example has been included in order to show that, even if 
we run into operational expressions of unfamiliar kinds, they will 
probably yield to the famihar rules. In 
some cases, the expert is able to detect y 
that the series which looks so fearsome v 
corresponds with some well-known and 
oft-tabulated function; thus, in the 
present example, we may write 



or else 


V ~ 2F 


0 i 


dt, 


where are functions called ‘Bessel Functions of 

orders 1 and 2 whose numerical values may be found in mathe- 
matical tables. The functions w^hich correspond with a few" of 
these more recondite ^-expressions have been included, for 
completeness, in Appendix I; the reader may obtain the second 
of the above expressions for v, by the use of forms (1) and* (16). 

7*7. With this example, the reader has been conducted to the 
very edge of our allotted territory. Like Alice on the mantelpiece, 
we stand with one foot in Looking-glass Land, where p behaves 
in strange ways and needs treating strangely. Alice got into 
trouble with Looking-glass arithmetic — ‘ “ She ca’n’t do sums a 
hitV the Queens said together, with great emphasis’ — and we 
might be no better if we ventured further. This book is not a 
guide to Looking-glass Land; in plainer words, it is not a guide 
to the behaviour of all circuits, but only to those simple, useful 
circuits wMch are mentioned in § 1*1. But I hope that the reader 
has found that even this restricted territory is varied enough to 
be interesting, for there is no doubt at all about its usefulness. 
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APPENDIX I 


OPERATORS— STANDARD FORMS 


The inclusion or exclusion of forms in this list must needs be 
somewhat arbitrary. The majority of useful forms are interpreta- 
tions of expressions like F{^) . 1 (or F{p ) . 1 ) ; and for the backbone 
of the list, F{p) covers a more or less complete series of the partial- 
fraction forms given in Chapter iii. A few derived forms clamour 
to be introduced, because of the simphcity of their ^expressions; 
these have their reference numbers inset, and distinguished by 
the letter {a). To interpret his partial fractions, the reader should 
look through the main sequence of numbered forms, until he 
finds those which he requires; but he should glance at the 
neighbouring (a) forms, in case his problem is one which termi- 
nates in this happy simplicity. 

The following have been excluded: 

(1) Derived forms which do not lead to simple ^-expressions; 
thus, from 






r . 1 = cos (Jl)t 


1 


;.l = 


1 — cos 0 )t 


Ct) 


,2 


we may deduce . o 

but we do not list this, since is not a recommended 

partial fraction. 

(2) Forms designed to enable lazy people to avoid making 
partial fractions, like 

:.l. 


p 


{p + a)(p + fiy 

If the reader wants an extended list with these forms added, 
he can easily compile one for himself. 

The operational equivalents of certain Besselfunctions, useful in 
examples like the one in § 7- 6, have been added. There are also a few 
operations of the form ^ 

m, 


p + a' 


to facihtate the use (where convenient) of the method given in 
§ 2*7. Many forms have been left out altogether, because they lie 
right outside the scope and methods of this book. 


OPEE.ATOES — STAWDAED POEMS 

Basic Poems 


(1) 

•w> 

II 

PH 1*^ 

(2) 

P-fit)=f'{t), if/(0) = 0. 


Opeeations epon 1 (OE 1) 

(3) 

1.1 =. 

p 

(4) 

^\-i 

II 

pH 

(5) 

■ ^ . 1 = e-“'. 
p + a 

(6) 

^ 1- e-«^ 

(j3 + a)»‘ («-!)! ■ 

(7) 

p2 

« . « . 1 — COS (Ot. 

(8) 

p 1 . 

(9) 


(10) 

p 1 1 , . 

(p^ + 0^2)2 • ^ ^ ~ COS 0)t). 


(10a) 

( 11 ) 

( 12 ) 

(12 a) 


^3 _ 


. 1 = < cos (j)t. 


.1 


-od 


P 

{p + ocf + w^ 


( COS (wi — — sin ojtl. 
0) 


p^ + ap 

— — - 1 = e-‘^^GOBO)t 
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7-2 


100 


(13) 
(13a) 

(14) 
(14a) 


OPBEATOES— STANDARD POEMS 

p + a / 


a/( 25^ 


f 




. 1 = Jo(a^). 


- i? (p-^J(p^-a^)Y' ■. _ r 

^(p^-a^)[. a ;• 


{at). 


,, ■ 1 = ^(aO- 

A^{p^-a?) 


(15) 1 = 

'JiP +^°‘-P)\ * 




(15a) 

(16) 

(17) 


P 


- 1 = e-“'/o(aO. 


+ 2ai)) ■ 


l p-^l(p^-cc^) y 


(18) 

(19) 

( 20 ) 
( 21 ) 
( 22 ) 

(23) 

(24) 


p + 


Operations of form 


-./(«) = e-‘‘f{\^>fit)dt. 


1 = j4(a<), 


■fit) 


23 + a 


{p + a} 


L_ . = e~“^ J . . . (w times) e“(/(if) dtdt .. 


1 1 

.1 = 

+ a a 

1 ;.« = iM-(l-e-0]. 


;p-ha a^ 
- . 6““^' 


^ + a 

1 


;P + a 
1 

p + a 


.cos = 

sinwi = 


a — /? 

a cos 6>^ + C£> sin cot — ae~^^ 
a sin (ot — (o cos + C£>e“®^ 


?i> 0. 

0 . 
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Forms (18) to (24) may be regarded as special cases of Heavi- 
side’s Shifting Theorem, namely 

(25) Fip + a).f{t) = e-‘^F{p).{e-tf{t)}. 

Note. In forms (13) to (17), J^, and are tlie standard 
notations for Bessel functions, which, may be found tabulated, 
for example, in Dale’s Five-Figure Tables of Mathematical 
Functions. 

Observe that (1) and (2) enable us to deal with forms which are 
standard except for a p more or less in the top or bottom of the 
j9-expression. Thus, (1) in conjunction with (16) gives 
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THE BINOMIAL THEOREM 
The expansion of (1 in ascending powers of x is given by 




2 ! 


3! 

n(n— 1 ) ... (to — m-l- 1 ) 

-| Qifn -j- , ^ ^ ^ 

m ! 


Note that n may be negative or fractional — ^in fact, in operational 
work, it almost certainly will be. The number of terms in the 
series will then be infinite. 

When a; is a number, the statement of the series must be hedged 
with provisos about the permissible values of x. In our work, 
however, x will be an operator, such as a/p; and the validity of 
the method must then be judged by results. However, if the 
p-expression tends to a finite value when^ is treated as a number 
and is made to tend towards infinity, binomial expansion of the 
j9-expression in powers of 1 jp will lead to a ^-series which means 
something. 
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TRIGONOMETRICAL AND OTHER FORMULAE 

This Appendix and the succeeding one are intended to provide 
a list of those formulae which are most likely to be useful in 
operational circuit work, not a complete list for all occasions. 
The hyperbolic functions (cosh, sinh) are therefore very scantily 
represented, because they can always be replaced by exponentials 
(g±aa;)^ and the latter have been preferred in the text of the book. 


(1) 

/y*3 /v*ti 

„ - w iV vV 

= ‘ + * + 11 +3! + - + i! ■"■■■■ 

(2) 

£;}x _j_ ^-jx ^2 ^4 

cos.i= 2 =1 2! + 4r 


(3) 

^ix __ ^-jx ^3 ^5 

sma;= =x-- + - 

x'^ 

(4) 

cosho; = ;r = 1 +-:r-. +“1 + 

2 2! 4! 

x^ 

(5) 

. , e® — e~® a;® 

smhx= 2 =* + ^1 + 5! + 

x^ 

^! + - 

(6) 

GOB^X-\-Bm.^X = 1. 


(7) 

cos 2a; = cos^o^ — sin^o; = 2cos^a; — 1 


= 1 — 2sin®a:. 


(8) 

sm2a: = 2sina;cosa;. 


(9) 

Qos{x+y) — cos a; cos 2 / — sin a: sin t/. 



cos (x - 1/) = cos X cos y + sin x sin 


(10) 

sin (x + 2/) = sin X cos 2/ + cos X sin 
sin(x-y) = sin x cos 2 /— cos x sin 2 /. 


(11) 

a cos X + 6 sin X = + 6®) cos (x - 



where cosji = a/^l(a^ + b^), 
acosoj + ftsina; = + sin(jr + ^), 

where cos f = b/^{a^ + b^), sin f = a/^{a^ + b^). 
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(12) cosa;cos2/ = J-[cos(a;-2/) + cos(a; + ?/)]. 

(13) sinccsin?/ = |[cos(a;-2/)-cos(a; + ?/)]. 

(14) sinajcosj/ = |-[sin(a;-j/) + sm(a; + «/)]. 

(15) cosaj + cosy = 2cos^-^cos^^^. 


(16) cosa;- cosy = -2sin^^ sin ^ 

2 2 

(17) sinaj + siny = 2sin^^cos^^. 

2 2 

(18) sina;-smy = 2cos^^sin^^. 

2 . 2 


The following Bessel function formulae may perhaps be 
inserted here: 


(19) 






(x/2p 


■ + ■; 


(x/2)^ 


(n+1) 2l(n-hl)(n-i-2) 

(xl2) 


3 1 (ti + 1 ) (ti + 2) (fi + 3) 


+ 




: + 


(x/2)^ 


(n-hl) 21 (n-hl)(n+ 2) 


+ 


(x/2f 


3 1 {ti -j- 1 ) (??' 4" 2) (ti + 3 ) 


+ 




2n 


(21) ^ Jni^) = 4-l(*) + Jn+xi^)- 


2'i} 

(22) -4(^)=4-x(x)-4h-i(^). 
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DERIVATIVES AND INTEGRALS 

These lists are severely limited to functions which, are useful in 
operational circuit work. The independent variable has been 
called t. 

Table oe Derivatives 


Function 

( 1 ) MV (product of two func- 
tions of t). 

(2) f(M), where u = (j){t) 
(function of a function, 
like sin ( 0 )^ 4 - 5)). 

(3) {at)\ 

(4) 

(5) co^o)t. 

(6) mio)U 

(7) cosha^. 

(8) sinha^. 

(9) e®^cos 
e“^sin 


(10) JJ^od). 

(11) Jo{<a). 


Derivative 
dv du 

ae^K 

— sin cot, 

0) COB (i)t. 

a sinh at, 
a cosh at, 

Aj{a^ + o)^) e^^cos {o)t + d + <p). 
^{a^ + 0 )^) sin (wt + d+ip), 

where cos §5 = a/^J{a^ + oj^), 
sin^ = oj/ljla^ + ci)^), 

n 

aJ„_i(a«)-- J„(ai). 

Tlf 

-J„(ai)-aJ„+i(ai). 

— aJi{at), 
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(12) 4(ai). 


c^In-x{ai)-~In{od,). 

Yh 


(13) Meet). 


Table 


Function 

dv 

(1) (product of 

two functions of t, 
one of which can 
be integrated). 

(2) { atr . 


ocIj{at) 

OF Integrals 

Integral (Indefinite) 

uv— Iv—dt, 

J dt 

(Integration by Parts.) 


O^npi^l 

n-Yl 


(3) 


(4) cos (i)t, 

(5) sin (i)t. 

(6) cosh at. 


(7) sinha^. 


—sin o)t. 

CO 

1 

cos (Ot. 

(0 

-sinha^. 

a 

- cosh at. 


(8) cos (cot + 0). 


sin (cot -i-0). 


i [a cos (o^ + <9) + a> sin (cot + ^)], 


^(a^H-o^) 


cos (cot + ff — ^). 


a ^- h(o 


z [a sin (cot + 6)-(0 cos (cot + 0)], 


where cos^ = al^(a^ + o)^), 
sin^ = (i)^{a^ + o)% 
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Function 

(9) OOS(t)jtOOS(Jt)2t. 
sin o)-^t sin (L> 2 t. 
sin o)itco^ 


(10) 


in-i • 


( 11 ) 


4(«0 

in-X * 


(12) 4(a«). 


Integral {Indefinite) 

Use trigonometrical formulae num- 
bered (12), (13), (14). 


4-iM) 
ocP-i ■ 

4-i(<^^) 

2 

-{4+i(a«) + 4+3(a«) 

+ 4 + 5(<^0 + • ■ ■ 

2 

-{4+i(a0-4+3(a^) 

+4+5(a«)- ...ad inf.}. 


(13) Uai). 
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THE SOLUTION AND PROPERTIES OP 
ALGEBRAIC EQUATIONS 

This is a bare summary of the properties of algebraic equations — 
just enough to enable the reader to find the factors of his 
^-expressions. For proofs of the assertions which are here 
dogmatically made, he should consult some such book as 
F. W. TurnbulFs Theory of Equations. 

Defikitioxs 
An expression of the form 

ax'^ -f + cx'^~^ + . . . + Za; -}- m, 

which (for brevity) we shall call F{x)^ is said to be ^polynomial of 
degree ninx. a,b, ...,m are constants, known as the coefficients of 
the various powers of x. We shall assume that these coefficients 
are real, and that a is positive. 

Any value of which makes F{x) = 0 is called a root of the 
equation F{x) ~ 0, 


Peopebties of Equations 

(1) An equation of degree n has n roots, which may be real or 
complex. Two or more roots may be the same. 

(2) When the roots x-^, ..., x,^ are known, F{x) may be 

written as the product of n factors, thus: 

ax '^ ... = a{x — x-f) {x — xf ) ... {x — x^). 


(3) This equation demonstrates the relations which exist be- 


tween the coefficients a, 6, 
namely 


m and the roots x- 


l3 •^2? 


bja = {x^ + (^2 + . . . x^) 

= - (Sum of roots). 


cja = +{x^X 2 + XjX^+ ...+X 2 Xq-{’ ...) 

= -f (Sum of all possible products of pairs of roots). 


mja = {-'l)^XjX 2 ^.^x^ 

= ( — 1)^ X (Product of roots). 
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(4) It follows that if we change alternate signs in F{x), thus: 

ax'^ — bx'^~^ -f cx'^'^^ — dx'^'~^ + ... + (— 1 = 0, 

we obtain an equation whose roots are numerically the same as 
the roots of F{x) = 0, but with their signs changed. 

(5) It also follows that if we multiply the coefficients by 

1, A, thus: 

ax'^ + bXx'^~‘^ -f- + . . . + mA^^ = 0, 


we obtain an equation whose roots are Xx^, Xx^, A^* 3 , etc. 

(6) Since the coefficients are assumed to be real, the complex 
roots (if any) can only occur in conjugate pairs, /? ± jy. 

(7) Hence any equation of odd degree must have at least one 
real root. 


(8) Hence also F{x) can alwa 3 ^s be split into real factors, of 
degree not greater than 2. 

(9) If F{x) is positive when x is replaced by one number g, but 
negative when x is put equal to another number r, then the 
equation F{x) = 0 has at least one real root between q and r . 

The last theorem affords a means of quickly finding the ap- 
proximate values of the real roots of an equation. For example, 

consider o ^ ^ 

x^-2x-5 = 0 , 


an equation which has become traditional as an instance in the 
theory of equations, ever since it was first propounded by John 
Wallis (1616-*1703), of Oxford. Let us jot down the values of 
F{x) for a few simple values of x : 

X — CO — 3 — 2 — 1 0 +1 “1-2 -h 3 +00 

F(x) -00 -26 -9 -4 -5 -6 -1 +16 +oo 

It will be seen that there is at least one root between 2 and 3, 
seemingly much nearer to 2. Actually this equation has only one 
real root, viz. 2*09455 .... 

Hoeneb’s Method oe obtaining Real Roots 

Having got somewhere near the root in this way, we can get as 
many decimal places as we require by the following method; 
some people are afraid of it, but its terrors are much exaggerated. 
It win be explained with reference to a cubic equation, but it is 
applicable to equations of any degree. 
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Suppose that inspection has shown the equation 

ax^ + bx^-\-cx-\‘d ^ 0 

to possess a root between 1 and 10, whose first digit is n. The 
coefficients are first written with their appropriate signs in a 
horizontal row, followed by the first digit of the root: 

® ^ 0 d {% 

The first coefficient a is now multiplied by n, and the product 
added to 6; this sum is multiplied by and the product added 


to c; and so on, thus: 




a b 

c 

d 

( » 

na 

nA 

nB 

(Sum) A 

{SumjE 

(Sum) G 



The process is repeated, stopping short at the last stage but one; 
then again, stopping short at the last stage but two; and so on, 
until the last process is but a single addition of na. The calculation 
for the cubic equation now looks like this : 


b 

c 

d 

na 

nA 

nB 

1 

~B 

C 

na 

nD 


'd 

E 


ma 




F 


The lowest digit in each column is usually brought into promi- 
nence by a dividing Mne ruled as shown; these constitute a new 
set of coefficients, namely a, F, F, C. 

It can be shown that the equation 

ax^ + Fx^ + Ex + C — 0 

has roots which are smaller by n than the corresponding roots of 
the original equation; so that a root of the original equation 
having the first digit has a counterpart root in the new equation 
which is less than 1. ' . 
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Tke equation 

ax^ + lOFx^ + lOOEx + lOOOC = 0 

has all its roots 10 times as large as the corresponding roots of 

' ax^ + Fx^ + Ex+ C — 0 (Theorem (5) of the last paragraph). ’ 

If we multiply the coefficients of the derived equation by 1, 10, 
100, in this way, we obtain an equation in which the root we 
are chasing (or its counterpart) again hes between 1 and 10. The 
whole process can be repeated, to obtain the second digit of the 
root of the original equation. 

The following numerical examples show how the work should 
be set out. 

Examples of Horner’s Method 

Example 1. To find that root of — 5a;^+l = 0, which lies 
between 2 and 3. 

2 ~5 

4 

4 

4 

70 

8 

8 

"lb 
8 

940 
2 

942 

2 

944 
2 

946 

From this stage, several more figures may be quickly obtained 
by a process analogous to the ' contracted division’ of elementary 


0 

-2 

6 


400 

312 

712 

344 


105600 

942 

106542 

944 


107486 


1 

-4 


( 2*41 


■3000 

2848 


152000 

106542 


-- 45458 
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arithmetic. Instead of adding O’s to the coefficients, we strike 
off one, two, three, . . . digits, but take account of the last digit 
neglected. The example therefore continues. as follows: 

//*2 107480 -45458 ( 2-414214 

37g 43146 


10786^ 

37g 


- 2312 
2165 


10824 

- 147 

2 

108 

10820 

39 

2 

43 

1023 


This result is correct to six decimal places. 

The degree of accuracy here obtained is greater than is usually 
needful in engineering problems. It has been thought well to go 
to this accuracy, in order to present a thorough illustration; but 
much less arithmetic is involved if the contracted method is 
begun after two stages, and the result is only wrong by 1 part in 
24,000: 


— 0 
4 

4 


0 

-2 

6 


1 

-4 


2-4143 


•3000 

2848 


11 % 


3 

4 

400 

312 1 

- 152 
106 

70 

712 


- 46 

8 

344 ' 


43 

78 

1050 

- 3 

8 

0 

3 

86 

1060 


8 

0 



IIH 


m 


Example 2. In § 7-4, we found it necessary to solve the equation 
216^c^— 162a;^— 60a: + 47 = 0. 

From the way in which the equation was derived, it appeared that 
we were interested only in a real root greater than |f, or 0-694. 
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The usual method of root-location shows that the equation has 
one negative and two positive roots ; both positive roots are less 
than 1, but one is a little greater than 0-7, while the other is less 
than 0‘6. We therefore write down the equation whose roots are 
those of the original equation multiplied by 10, namely 

216a:3-1620a;2-6000x + 47000 = 0, 


and seek that root whose first digit is 7. 

The Horner calculation, using the contracted method after two 
stages, looks hke this (the coefi&cients have all been divided by 4) : 


lU 


-405 -1500 11750 ( 0-7087 

378 - 189 -11823 


- 27 -1689 -73000 

378 2457 66131 


351 76800 

378 586^ 


- 6869 
6236 


7200 8266^ 

4 : 5890 


633 


733 


8850 

52 


737 


8900 


u 


This degree of accuracy is sufficient. If the reader follows the 
working through for himself, he wiU easily discover what happens 
at the stage when the digit is zero. 


The Mechanism of Hoenee’s Method 

The process of solution wiU be much better understood if Ex- 
ample 1 is further considered. 

The roots of the equation are actuaUy 

-0-5858, +0-5000, +2-4142, 

to four places of decimals. 

The first reduction makes them 

-2-5858, -1-5000, +0-4142. 

(Notice that the product has changed sign. This is to be seen in 
the last coefficient of the equations.) 
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The first imiltiplication gives 

-25-858, -15-000, +4-142. 

The second reduction gives 

-29-858, - 19 . 000 , +0-142. 

The second multiplication gives 

— 298-58, — 190-00, + 1*42, and so on. 

It is appaient that the root which we are chasing is rapidly 
becoming by far the smallest root of the equation, and the reduc- 
tion process cannot cause any of the roots to change their signs. 
This gives us a rule for determining the successive digits (after the 
first, which was found by location of the roots): Find the largest 
digit which does not cause the last coefficient to change sign. 
(In exceptional cases, where there are two roots very close to- 
gether, this may break down; the reader can try out hypothetical 
cases to test this. The root-locating process may then have to be 
repeated.) 

To find a negative root, change the signs of alternate coefficients, 
so that the root becomes positive. 

To find a root not lying between 1 and 10, multiply or divide 
all roots by such a multiple of 10 as will bring it within those 
limits. 

The Quadratic Equation 
The roots of the quadratic equation 

+ c = 0 

may be found immediately from the formula 


^2 ~ 


— — 4ac) 

2a 


which is given in elementary algebra text-books. 


The Cubic Equation 

Every cubic equation must have at least one real root. This may 
be found by Horner’s method, and the equation reduced to a 
quadratic by dividing by the corresponding factor. 


GET 


8 
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The Quabtic (Biquadratic) Equation 

A qiiartic equation may have four real roots^ or two, or none. 
The last ease presents us with a new difficulty, which may be 
removed in the following manner. 

We know that it must be possible, in at least one way, to re- 
solve an expression of form 


ax^ -r bx^ + cx^ + dx-\-e 


into real quadratic factors. Therefore it must be possible to find 
real numbers/, g, /^, such that 

(Note that the coefficient of re in the first squared expression has 
been assigned that value which leads, on squaring, to the right 
coefficient for x^,) 

Equating coefficients, we find 



c 

a 





d 

a 



e 

a 


P-hK 


The first equation gives 




the third, = 

a 

Substituting in the second, 


which is a cubic equation in/, and therefore within our compass. 
An example of this method will be found worked out in § 7-4. 


115 


ALGEBRAIC E QU ATIO'JTS ■ 

The Qhintic Equation 

Every qiiiiitic equation must have at least one real. root, and can 
therefore be reduced to a quartic. 

Equations of Higher Order 

The sextic equation' may have three pairs of complex roots. If it 
has, the methods so far given all break down. Methods are known 
for overcoming this obstacle, but an engineer is unlikely to find 
himself repaid for the labour of learning or of applying them; 
therefore I deliberately refrain from carrying this discussion 
further. Interested readers will find the usual method described 
in Heaviside' s Operational Calculus, by E. J. Berg; it is called 
Graeife’s method. 

Conditions for Real or Complex Roots 

(а) Quadratic, ax^^bx-\-c = 0. 

If ~ 4ac > 0, 

the roots are real; if this expression is less than 0, they are com- 
plex. 

(б ) Cubic , ax^ + bx"^ + ca: + c? ~ 0 . 

If 21 — bh^ + 4- 4d6® — 1 Sabcd < 0, 

aH the roots are real; if this expression is greater than 0, there is 
one real root and a pair of comxilex roots. 

(c) Biquadratic, ax^ -{■bx^ + cx'^-i-dx + e = 0. 

If 4:{l2ae—3bd + c^f 

— {12ace + dbcd ~ 21 ad^ — 21 eb^ — 2c^)^ < 0, 

the equation has two roots real and two complex; if this expres- 
sion is greater than 0, the roots are all complex unless both the 
following conditions are also satisfied: 

8aC”36^<0, 

16a2(12ae 4*0^) - (8ac -- < 0; 
in which case the roots are all real. 
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CONDITIOKS FOR STABILITY 

The following conditions must all be satisfied if the equations are 
to have no roots with a positive real part. Applied to the bottom 
line of a j^-operator, these conditions are sufficient to ensure that 
terms like or e^^coscot shall not occur. They are given by 
E. J. Routh, Advanced Rigid Dynamics, pp. 192-202 (1892 
edition), and are usually known as ' Routh’s Criteria of Stability 
It is assumed that the first coefficient a is positive. 

(а) Quadratic, + + c = 0. 

Conditions for stability are 

6>0, oO. 

(б) Cubic, ax^ + bx^ + + cZ = 0. 

Conditions for stability are 

6>0, bc — ad>0, cZ>0. 

(c) BiquadjTatic, ax^ + bx^ + cx^ + dx + e = 0. 

Conditions for stability are 

6>0, bc — ad>0, bcd — ad^ — eb^>0, e>0. 

Conditions for Existence of Complex 
Roots OF FORM -0(1 ±j) 

These are the conditions discussed at the end of §6-8; applied 
to the bottom line of a j5-operator, they foreshadow the presence 
of oscillatory terms like cos wt or sin o)t, 

{a) Quadratic, ax^-\-bx-\’C — 0. 

The condition is 

6^ — 2ac==0. 

{b) Cubic, ax^-^’bx'^ + cx + d ^ 0. 

The condition is 


6V+2ac^-f 2<I6^ — 4a5ctZ = 0. 
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